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SUMMARY 

The  Brownian  motion  of  a  test  body  due  to  multiple  interactions  with 
field  particles  is  investigated  within  or  almost  within  the  framework  of  Markoff's 
processes.  First,  Markoff's  processes  are  studied  as  presenting  such  multiple 
interactions.  Based  on  the  study  and  by  means  of  Markoff's  method  of  random 
flights,  we  investigate  the  Brownian  motion  of  an  elastic  test  body  submerged  in 
a  rarefied  gas  constituted  of  elastic  molecules,  under  the  condition  that  mutual 
interactions  among  field  particles  are  negligible.  It  is  shown  that  there  is  no 
difference  in  effect  between  temporal  repetitions  of  random  binary  collisions  and 
multiple  collisions  (random  binary  collisions  superposed  at  one  moment  of  time), 
so  far  as  the  friction  and  diffusion  of  the  test  body  in  momentum  space  are  con¬ 
cerned.  The  situation  is  similar  when  a  test  body  with  electric  charge  is  sub¬ 
merged  in  an  electron  gas,  i£  the  mutual  interactions  among  electrons  are  ignored. 


h  This  research  was  supported  by  the  Office  of  Naval  Research  under 
Contract  No.  Nonr  839(38),  Project  No.  NR  0t>l-135. 

*  Visiting  Professor,  Department  of  Aerospace  Engineering  and  Applied 
Mechanics. 


\ 


It  is  not  feasible,  however,  to  ignore  those  mutual  interactions  of  field  electrons 
and  to. represent  electronic  multiple  interactions  by  temporal  repetitions  of 
random  binary  interactions,  each  of  which  takes  place  independently:  Fluctua¬ 
tions  of  limitlessly  large  amplitudes  in  the  spatial  distribution  of  electrons, 
which  may  possibly  take  place  in  this  approximation, do  not  seem  realistic,  be¬ 
cause  a  limitless  concentration  of  potential  energy  accompanying  a  concentration 
of  electrons  in  a  local  spot  cannot  be  permitted.  Amplitudes  of  such  fluctuations 
and/or  microscopic  disturbances  must  have  a  certain  maximum  limit.  [The 
situation  does  not  change  even  when  the  interaction  force  law  is  of  the  Debye- 
Huckel  type.  ]  A  kinetic  theoretical  scheme  of  treating  fully  ionized  gas  in  the 
light  of  this  fact  is  proposed. 


u 


TABLE  OF  CONTENTS 


Section  ^*8* 

I  Introduction .  1 

II  Collisions  as  Markoff's  Processes .  4 

III  Elastic  Test  Bodies  in  a  Rarefied  Gase .  13 

IV  An  Electrically  Charged  Test  Body  in  an  Electron 

Gas .  24 

V  A  Scheme  of  Kinetic  Theoretical  Treatment  of  an 

Electron  Gas .  31 

VI  Concluding  Remarks .  35 

VII  References .  36 

Appendix  A.  The  Derivation  of  the  Fokker- Planck 

Equation  from  the  Smoluchowski  Equation  .  37 

Appendix  B.  The  Fokker-Planck  Equation  Derived 

from  the  Boltzmann  Equation .  40 

Appendix  C.  A  Solid  and  Elastic  Test  Body  in  a 
Rarefied  Gas  Constituted  of  Solid  and 

Elastic  Molecules . 44 

Appendix  D.  Markoff's  Method  of  Random 

Flights .  57 

Appendix  E.  Multiple  Collisions  of  an  Elastic 

Test  Body  in  a  Rarefied  Gas .  58 

Appendix  F.  Useful  Integrals 

iii 


71 


SECTION  I 


INTRODUCTION 

The  difficulty  in  kinetic  theoretical  treatments  of  an  electron  gas  or 
a  fully  ionized  gas  may  be  attributed  mainly  to  the  fact  that  we  are  not  able  to 
treat  precisely  more  than  two  body  problems.  In  Part  I  of  this  report 
we  showed  that  methods  of  treating  multiple  interactions  among  electrons  in 
an  electron  gas  in  accordance  with  the  B-B-G-K-  Y  hierarchy  are  not  plausible. 

In  the  present  part,  we  attempt,  by  means  of  several  "mental  experiments", 
to  find  characteristic  effects  which  distinguish  multiple  interactions  from  superposed 
binary  interactions.  Let  us  first  suppose  that  the  effect  of  multiple  collisions 
between  a  test  body  and  field  particles  is  divided  into  three  parts:  (i)  the 
effect  of  the  forces  exerted  on  the  test  body  simultaneously  by  many  field  par¬ 
ticles,  (ii)  the  effect  of  reactive  forces  exerted  by  the  test  body  on  field  par¬ 
ticles,  (iii)  the  effect  of  mutual  interactions  among  field  particles.  Since  we 
do  not  know,  in  general,  the  way  to  synthesize  the  total  effect  from  those  par¬ 
tial  effects,  that  is,  the  many-body  problem,  the  classification  seems  simply 
conceptual  with  no  definite  physical  meaning.  We  note,  however,  that  it  is 
possible  to  concrive  particular  conditions  under  which  effects  (ii)  and  (iii)  are 
ignored  either  completely  or  partly,  and  yet  interactions  are  multiple.  By 
studying  such  particular  cases,  we  may  obtain  some  aspects  of  statistical 
effects  of  the  multiplicity  of  interaction,  and, we  hope, may  find  some  clue  for 


considering  multiple  interactions  among  charged  particles  in  a  fully  ionized  gas. 

We  suppose,  first,  that  a  test  body  is  an  elastic  body  with  linear 
dimension  D  which  is  sufficiently  larger  than  the  average  distance  of  two 
neighboring  molecules  of  the  gas  in  which  the  test  body  is  submerged.  Further, 
we  suppose  that  D  is  much  shorter  than  the  mean  free  path  of  molecules  of  the  gas, 
and  that  the  mass  of  the  test  body  denoted  by  M  is  much  larger  than  the  mass 
of  a  molecule  m.  Under  these  conditions,  the  interactions  between  the  test 
body  and  field  molecules  may  be  multiple.  However,  the  mutual  interactions 
among  field  particles  are  negligible.  We  investigate  the  statistical  behavior 
of  the  test  body  by  assuming  that  the  gas  is  in  thermal  equilibrium.  At  first 
glance,  it  seems  possible  to  treat  the  particular  problem  by  means  of  the  usual 
technics  of  kinetic  theory  based  on  the  Liouville  equation  of  those  particles. 

We  note,  however,  that  collisions  are  strong  regarding  field  molecules  even 
though  the  collisions  are  weak  for  the  test  body.  Under  this  circumstance  the 
coarse-graining  of  the  Liouville  equation  deriving  the  statistical  behavior  of  the 
test  body  is  still  complicated.  We  wish  to  find  some  other  approach; 

We  see  that  the  collision  processes  under  consideration  may  be  Markoffian  or 
almost  Markoffian.  Markoff  processes  have  been  well  studied  as  a  mathematical 
theory.  *  The  Brownian  motion  of  a  test  body  in  a  fluid  is  well-known  as  a 
typical  example  of  the  theory.  Usually,  however,  the  field  particles  are  not 
objects  of  direct  investigation  in  the  theory:  Physical  interpretations  of  the 

-  ^ 

Since  Einstein's  pioneering  study  of  the  Brownian  motion  in  the  beginning  of 

this  century,  the  theory  has  been  developed  by  many  authors.  The  names 

of  Smoluchowski,  Fokker,  Planck,  Ornstein,  Burgers, Fiirth,  Uhlenbeck,  and 

Chandrasekhar  are  well  remembered.  2,  3,  4,  5 
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results  are  achieved  by  providing  data  of  relevant  physical  quantities  such  as 

2  3 

friction  and  diffusion  tensors  from  outside  the  theory  ’  .  For  instance, 

the  friction  coefficient  is  often  provided  by  Stokes'  friction  law.  One  of  the 
primary  interests  in  those  studies  has  been  to  show  that  the  distribution  of 
the  state  of  a  test  body  finally  becomes  Maxwellian.  It  is  noted,  however, 
that  the  theory  itself  has  not  proved  the  assertion  that  the  temperature  of 
the  final  state  of  the  test  body  is  the  same  as  that  of  the  fluid.  * 

In  this  study  by  taking  advantage  of  the  simplicity  of  collision 
mode  as  stated  above,  we  attempt  to  include  the  dynamics  of  those  collisions 
within  the  framework  of  the  theory  of  Markoff  processes  as  stated  in  Section  II. 
In  Section  III  we  first  consider  a  special  case  where  the  test  body  and  the  field 
molecules  are  extremely  rigid  so  that  the  period  of  each  collision  is  extremely 
short  and  hence  collisions  are  binary.  It  is  shown  that  the  temperature  of  the 
final  distribution  of  the  test  body  is  the  same  as  that  of  the  gas.  **  We  then 
assume  that  the  period  of  each  collision  is  finite  so  that  many  collisions  occur 
simultaneously.  Physical  quantities  of  interest  are  obtained  by  means  of 
Markoff's  method  of  random  flights  under  various  conditions.  In  view  of  the 
conclusion  attained  in  this  section,  we  consider  some  aspects  of  multiple 

❖ 

According  to  the  principle  of  statistical  mechanics,  we  may  say  so.  The 
question  is  open,  however,  if  the  test  body  is  much  larger  than  molecules. 

## 

Related  to  the  investigation  in  this  section,  it  is  shown  that  the  Boltzmann 
equation  which  is  obtained  under  the  assumption  that  collisions  are  binary 
and  Markoffian,  is  reduced  directly  to  the  Fokker-Planck  equation  by  assuming 
that  interactions  are  weak.  See  Appendix  B. 
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interactions  in  an  electron  gas  in  the  following  section.  Finally,  a  set  of 
basic  equations  for  electron  gas  is  proposed.  The  detailed  analyses  of  the 
equations  will  be  given  in  Part  III. 


SECTION  II 

COLLISIONS  AS  MARKOFF'S  PROCESSES 
First  we  remember  Markoff's  proces  ses,  well  studied  in  the  mathema- 

-_4  —4 

tical  theory.  By  taking  p  for  a  variable  of  state,  t  for  time,  and  f^(p,  t)  for  the 

distribution,  we  assume  that  the  evolution  of  f  is  a  Markoff  process. 

o 

In  other. words,  f  is  assumed  to  satisfy  the  Smoluchowski  equation: 
o 

fQ(p,  t  +  At)  =  IJj  fQ(p',  t)  Kp.'t/p,  t  +  At)  d3p'  (2.  1) 


where  'll  is  defined  as  transition  probability.  By  putting 


p  =  p  -  Ap 


we  may  write 

'll  (p',  t/p,  t  +  At)  d3p'  =  -  4i(p  -  Ap,  t/p  -  Ap  +  A  p,  t  +  t)  d3  Ap 
=  cp  (p  -  Ap,  t;  p,  A  t)  d3  ,  p 
By  expanding  Cp  in  a  Taylor  series,  wo  obtain 


if  (p  t/p,  t  +  A  t)  d3  p  ' 

=  £  ¥  (P>  t;Ap,  At) 


.  "  S  4  ■  1  -  ,  - 

p-  df  +  17  p,  p:dJIp 


3! 


-  TT  L  PA  P-  P  ■  ^777 


1  dp  Ap  dp 


(2.2) 


(2.3) 


(2.4) 
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We  put  for  cp 

JJJcp(p,  t;  Ap,  At)  da  Ap  =  1  (2.5) 

— *  “ *  H 

If  we  define  <Ap>^,  <ApAp>^, - ,  by 

<Ap  >  =  J]Jap  cp  (  p,  t;  Ap,  At)  daAp, 

<  Ap  Ap>^t  =  JJJ Ap  Ap  Cpd3Ap,  (2.6) 


Eq.  (2.  1)  yields,  as  shown  in  Appendix  A, 


dfo(p,t) 

at 


At+  - 


_ o 

at8 


(At)8  + 


3! 


a 

ap 


•(<4pyo» 


i  a? 

2!  ap  ap 


a3  f 

^(At)3+  — 
(<A?A?>tfo) 


TT  Iplpf;:  (<4PAPApyo)  + 


It  is  noted  that  <Ap  >  ,  <ApAp>  - are,  in  general, 

<Ap>  /At,  <Ap  Ap>  /At,  .  .  .  have  certain  limit  values, 
At  At 

a  short  period  denoted  by  T,we  may  write  for  (2.  7) 

IT  =  '  a|' (<Ap  >fo)  *  2\  a^Tp:  (<^P^P>f0) 


IT  ayffTp  :  (^P  ‘P  -p>f0)  +  — ■ 


(2.7) 

functions  of  At.  If 

as  At  tends  to  decrease  to 


(2.8) 
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where 


<  Ap  >  =  lim 


<Ap  >At 


it  -T  At 


'>  (independent  of  T) 


<ApAp>  =  lim  Ap  At 


At  -T  At 


(independent  of  T) 


(2.9) 


If  they  fail  to  have  certain  limit  values,  it  is  difficult  to  give  proper  physical 
interpretations  of  Eq.  (2.  7). 

Secondly,  we  consider  the  same  process  as  a  physical  problem.  To 
this  end,  we  consider  a  test  body  whose  momentum  changes  as  a  function  of  time 
due  to  forces  exerted  at  random  from  the  outside.  By  considering  similar 
systems,  we  may  define  the  distribution  function  of  the  test  body  f^,  as  a 
function  of  its  momentum  and  time.  Each  of  the  external  foces,  denoted  by 
y^,  is  assumed  to  appear  with  its  proper  probability  w^  with  respect  to  a  unit 
time,  and  to  continue  for  its  proper  period  T  .  The  change  of  the  momentum 
of  the  test  body  due  to  is  given  by 

t  +  T 

Pr  =  Jr  r  yr  dt  (2.  10) 

t 

r 

where  t  is  the  initial  time  when  y  begins  to  appear.  We  assume  that  w^_  is 
uniform  during  an  elementary  period  At  which  is  longer  than 

At  >t  (2.  11) 
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We  also  assume 


w  t  <<1  . 
r  r 


(2.  12) 


In  this  section,  we  do  not  consider  the  physical  cause  of  each  force;  we  simply 
assume  that  y  and  w^  do  not  change  due  to  the  change  of  momentum  of  the  test 
body  during  At.  To  apply  the  result  of  the  present  section  to  physical 

problems,  it  is  necessary  to  choose  cases  where  the  assumption  is  satisfied. 

The  number  of  forces  which  occur  during  the  period  At  is  given  by 

N 

<V>At  =I  V*  (2*13) 

r  =  1 


where  N  is  the  total  number  of  forces  and  it  is  assumed  that  a  force  does  not 
appear  twice,  due  to  some  mechanism  of  the  external  cause.  It  is  possible  that 
a  force  exists  partly  within  the  period  At  and  partly  beyond  (before  or  after) 
the  period.  Due  to  (2.  11)  and(2.  12),  however,  the  number  of  such  forces  is 
negligibly  small,  otherwise  the  processes  due  to  these  forces  are  not  Markoffian. 
Under  those  conditions,  the  probable  value  of  Ap  ,  the  momentum  gain  of  the 
test  body  in  At,  is 


N 


<  Ap 


>  = 

At 


I 


p  w  At 


r  =  1 


(2.  14) 


It  is  noted  that  Ap^  fluctuates,  if  similar  experiments  are  repeated.  The 
probability  distribution  of  Ap^  is  denoted  by 


cp(Ap;  p) 

where  p  is  the  momentum  of  the  test  body. 


(2.  15) 
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Of  course  we  have 


HT  cp  (Ap;  p)  dJAp  =  1 


(2.  16) 


According  to  the  definitions  of  f  and  cp,  we  write 

o 

fQ(p.  t  +  At)  =  JJJ CP  (  Ap;  p-Ap)  fQ  (p-Ap,  t)  d*Ap 


(2.  17) 


This  is  obviously  equivalent  to  Eq.  (2.  1),  and  hence  Eq.  (2.  17)  yields  Eq.  (2.  7). 

In  the  following,  we  consider  cp(Ap;p^  itself.  The  probability  of  no  force 
appearing  in  At  is 

N 

(2.  18) 


s(o)  =  n  (i  -  w  At). 

r=l  r 


The  probability  of  only  one  force  appearing  in  At  is 
'  N  N 

S(l)  =  S  s  =)  w  At  Fi  (l-w  ,At) 

4j  r  Lt  t  t,  r 

i  i  r*r 

r=l  r=l 


Similarly, 


r, 

rx  <  rs  rx<rs 


=  )  V  w  Atw  At  ,n  (l*w  At) 
L>  L  r  i  r2  r3/r1(  rs  r3 


ri  <r3  <-  -  -<r  n 

*  n 

= y  ^  —  y  w  w  —  w  (At)  r (i-w  At) 

L.  L-  L  rx  ra  r  r 


ri  <rz - <r 


(2.  19) 


rn+l  *  ri-r2'-"’rn 
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Case  I.  We  first  consider  the  case  where 


<v> 

At 


At  <  1, 


N  >>  1 

In  this  case,  we  may  expect  that 


(2.  20) 


s  =w  At  <<  1  (2.21) 

r  r 

and  the  higher  powers  of  w^  At  are  negligible.  Hence  S(0)  given  by  (2.  1$)  and 
S(l)  predominate.  We  may  put  for  <Ap  Ap>At 

<ApAp  >  =)  s  p  p  (2.22) 

At  L>  rrr 

Similarly, 

<AP  AP  Ap>'t  =  £  srPr  PrPr 

(2.23) 


In  order  to  satisfy  condition  (2.  21)  under  the  restriction  (2.  11),  it  is 
necessary  to  assume  that  either  ^  is  extremely  short  or  w^  is  very  small. 
When  a  rigid  test  body  is  submerged  in  a  rarefied  gas  composed  of  rigid 
molecules,  T  may  be  extremely  short  and  \\  may  be  small.  * 

*The  Boltzmann  equation  is  derived  under  the  same  conditions.  Hence  it  is 
expected  that  the  Boltzmann  equation  may  be  reduced  to  (2.  8).  This  is  shown 
in  Appendix  (B). 


1 1) 


Case  II.  If  we  assume 


<v>. 


At 


=1% 


At  >>  1 


(2.  24) 


we  have  to  consider  all  the  S's  given  by  (2.  18)  and  (2.  19).  In  this  case, 
<Ap>^  given  by  (2.  14)  is  still  valid:  we  have,  however,  for  <Ap  Ap>( 


At 


N 


<apsp>st  •  ^  .r  Pp  Pr  *11  «ri  ra  (pti  +  pta)  (prj  *  prj) 


*--*ll-i 

ri<r8<  <r 


ri  *  ra 


sr  r  ...  _  (P_  + - +Pr  )(p  +---+P  ) 

ri  fs - r  rx  r  rx  r 


(2.25) 


Considering  the  law  of  large  number,  it  is  expected  that,  as  <V>  increases, 
the  fluctuation  in  the  number  v  decreases*  and  we  may  simply  take 

<APAP>At  =  H'"I  %>...  r_(Pri+---+Pr_Mp'’ri+--*-+Pr  ) 

-  -<  r —  n  n  15 


n  ’  <V>At 


(2.26) 


We  may  consider  <ApAp  Ap>^  - in  a  similar  way.  Those  quantities  so  far 

obtained  do  not  necessarily  have  definite  physical  meanings:  For  those  quanti¬ 


fy  taking  At  extremely  long,  we  may  consider  the  forces  appearing  in  At  as  a 
"population".  In  the  present  argument,  we  are  doing  a  random  sampling  from 
the  population.  See  any  text  book  of  "Statistics". 


ties  to  have  certain  physical  meanings,  it  is  necessary  that 


<v> 


<v>  = 


At 


At 


,  <Ap>  = 


:Ap> 


At 


<Ap  Ap  Ap> 


At 


,  <  Ap  Ap  Ap>  = 


At 


At 


<Ap  Ap  Ap  Ap  >  = 


<ApApApAp> 


At, 


At 


(2.27) 


must  converge  to  certain  values  which  are  independent  of  At,  as  At  decreases 
toward  its  lower  limit  given  by  (2.11).  It  is  also  necessary  that 


at  this  limit. 


2?f 


df 


IF"  '  if  <<:l’ 


(2.  28) 


We  simply  assumed  that  w^  and  y  are  invariant  regardless  of  the 
sequence  of  appearance  of  more  than  one  force  in  At.  As  a  matter  of  physics, 
this  assumption  is  feasible  only  under  certain  conditions.  In  the  following 
sections,  we  will  pay  special  attention  to  this  assumption  in  each  case  so  that 
the  assumption  will  be  feasible. 
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SECTION  III 


ELASTIC  TEST  BODIES  IN  A  RAREFIED  GAS 

Our  first  question  is  as  follows:  Is  there  any  difference  between  the 
effect  of  binary  interaction  and  that  of  multiple  interaction  when  mutual  inter¬ 
actions  among  field  particles  are  ignorable?  In  order  to  answer  this  question, 
we  investigate  several  physically  conceivable  examples  within  the  framework 
of  the  restrictive  conditions  considered  in  the  preceding  section.  We  consider 
a  mechanically  elastic  test  body  of  linear  dimension  D  and  with  mass  M,  sub¬ 
merged  in  a  gas  composed  of  spherical  particles  (molecules)  of  one  species  with 
diameter  a  and  mass  m,  in  thermal  equilibrium.  We  assume 


M  > > m, 

(3.  1) 

1 

D  >>n  '  3  »a, 

(3.2) 

\ 

X>>  D 

(3.  3) 

where  X  is  the  mean  free  path  of  field  particles  (molecules),  n  the  number 

i 

density  of  field  particles.  It  is  noted  that  n  3  is  of  the  same  order  as  that 
of  the- average  distance  between  two  neighboring  particles.  By  (3.  1)  it  is 
assured  that  the  collision  between  the  test  body  and  a  field  particle  is  weak 
so  far  as  the  test  body  is  concerned:  by  (3.  2)  the  number  of  collisions  in 
'  t  may  be  large,  and  by  (3.  3)  the  effect  of  mutual  interactions  among  field 
particles  are  negligible.  Due  to  the  last  condition,  field  particles  are  molecular- 

1  1 


disorde  red. 


1.  Binary  Collision  (t  -0) 

Case  1  -a.  For  the  convenience  of  comparison,  we  first  consider 
binary  collision.  This  condition  may  be  realized  by  assuming  that  the  test 
body  and  field  particles  are  extremely  rigid  so  that  the  period  of  a  collision 
is  extremely  short  or  T^=0.  By  taking  a  sphere  of  radius  D  for  the  test  body, 
as  obtained  in  Appendix  C,  the  average  number  of  field  particles  colliding  with 
the  test  body  in  At,  which  satisfies  (2.  21)  ^w^  At  <1,  is  given  by 

&X)  „  ,3.4, 

At  m 


where 

R=  2  +  ° 

2  2 

n  =  number  density  of  molecules. 

The  average  momentum  given  by  those  particles  to  the  test  body  is 


(3.5) 


<AP>At  =  *  f  mv<v>/\t  (3,6) 

where  v  is  the  velocity  of  the  test  body  relative  to  the  gas.  Further  we 
obtain,  according  to  (2.  22), 


<Ap  Ap> 


At 


4  2.  2  k  T  ,  ,  .  r* 

’  3  m  '  —  l<V>it  f 


6  :  unit  tensor 

In  those  derivations,  it  is  assumed  that 
2 


mv 

kT 


<<1  . 


(3.7) 


(3.8) 
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This  assumption  may  be  justified  by  considering  (3.  1).  Similarly,  according 
to  (2.  23),  we  have 

I  Components  of  <  Ap  Ap  Ap  |  oc  ma  (  -jjj-  )  v  At, 

(3.9) 

(Components  of  <Ap  Ap  Ap  Ap^Joc  m4  (~~ )3 At 


By  considering  that 
df 


lf"o“0  l*r] 


d5* 

r=5T-  /f  -o  [T7T=r  +  (  r=.)8] 

dp  op  '  o  MkT  kT 


d3f 

.  /f  [  (  ^77-  +  )  (  r^r  )  ] 

dp  dp  dp  o  MkT  (kT)a  kT 


we  may  ignore  those  higher  order  moments  as  given  by  (3.  9).  According  to 
these  considerations,  Eq.  (2.  7)  yields 


4 

—  m  v 


l^!)f 

At  0 


1  da  .4  2  2kT 

2  dpdp  :  (  3  m  m 


2it 

At 


(3.10) 


It  is  easily  shown  that  the  right-hand  side  of  the  equation  vanishes  if  we  put 

M 


e  ,  ...  .3/2  r  M 

f  =  ao  (  2"nkT  )  exp  L'^ 


i! 


2kT  Ms 


]. 


(3.  11) 


M  v. 
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Case  1  -b^. 


For  the  convenience  of  later  comparison,  we  here 


take  for  the  test  body  a  rigid  cube  moving  with  velocity  v  which  is  perpendicular 
to  one  of  the  faces  of  the  cube.  In  this  case,  O  is  the  length  of  a  side.  The 
results  obtained  in  Appendix  C  by  taking  v=  (v,  0,0)  are 


„  ^  3  .  2kT .  1/2  „  A 

<v>a.=  7S<— >  ■> DS4*' 


(3.12) 


-  3mV<V>Af 


<ap  >  .  =  <ip  >  »  o 

y  At  z  At 


(3.  13) 


<Ap  Ap  >  =  <i)p  Ap  >  =  <Ap  Ap  > 

x  x  At  y  y  At  z  z  At 


4  _  .  2kT.  ,  ^ 

=  T  m  (  —  )  <v>Af 
3  m  At 


(3.  14) 


Since  Ap^,  Ap^,  and  Ap^  are  events  independent  of  each  other*,  it  is  natural 


th-t 


<  Ap  Ap  >  =  <Ap  >  <Ap  >  =0 

*x  *y  At  rx  At  *y  At 


(3.15) 


<Ap  Ap  >  =  0 

x  z  At 

Taking  advantage  of  the  simplicity  of  treatments,  we  obtain 


.  a  a  s  -32  3  .  2kT . 

<  Ap  Ap  Ap  >,  =  —  m  ( -  )  v  <  V  > 

rx  fx  At  3  m  At 


<Ap  Ap  Ap  Ap  > 
x  x  x  rx  At 


32  4  2kT  . 2  .  . 

—  m4  (  -  )3  <v> 

3  m  '  At 


(3.  16) 


*The  particles  which  contribute  to  Ap  are  different  from  the  particles  which 
contribute  to  Ap  . 

y 
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The  other  components  of  <Ap  Ap  Ap  >and  <Ap  Ap  Ap  Ap  >  are  shown  to  vanish.  By 
comparing  (3.  12),  (3.  13),  and  (3.  14)  respectively  with  (3.  4),  (3.  6),  and  (3.  7), 
we  may  say  that  the  above  two  cases  are  similar  to  each  other.  In  other  words, 
Eq.  (3.  10)  is  valid  in  this  case  too. 

Case  1-b..  On  each  of  the  faces  of  the  cube  which  are  parallel  to  v,  we 

assign  <v  /6  field  particles  (with  no  fluctuation)  as  colliding  on  the  face  in  A  t 

On  the  face  which  is  toward  the  direction  of  v ,  ^  +  JL  jya  n  ^t)  field 

particles  are  assigned*,  and  on  the  opposite  face,  /  — — —  .  —  o3n  At)  particles 

'6  2 

are  assigned.  The  distribution  in  the  configuration  space  is  then  no  longer  at 
random.  Assuming  that  those  particles  are  distributed  according  to  the  Maxwell 
function  in  the  momentum  space,  we  obtain,  at  the  limit  T^=At  =  0. 


<Ap  Ap  >  =  f  mkT<v> 

x  x  At  3  At 


(3.  17) 


<Ap  Ap  >  = - -  0. 

y  y  At 


Case  1-b^.  We  now  assume  that  the  distribution  of  field  particles  in  the 
configuration  space  is  at  random.  However,  each  particle  is  assumed  to  have  the 
same  magnitudes  of  momentum  components 


*See  Appendix  C,  These  are  the  average  numbers  of  particles  colliding  on  the 
two  faces  which  are  perpendicular  to  v  in  Case  1-b^, 
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(3.  18) 


j]Jm|  c^l  fd»c 


I  T  '*/2 

'  '  2n  ’ 


In  this  case  we  obtain  for  v=  (v,  o,  o) 


.  _  6nDa  ,  kj  .1/2 

<V>  =  — — r  12  \  ~  >  to  , 

At  (2n)  1/4  m 


<Ap  >  =  -  —  m  v  <v>  , 

*x  At  3  At 


(3.  19) 


<Ap  >  =  <  Ap  >  =0 

FyAt  Fz  At 


< Ap  Ap  >  =  <Ap  Ap  >  =  <Ap  Ap  > 
x  x  y  y  z  z 


I2i»J»s<v>m/3  :i-  mkT<n, 


Case  1-b^  •  Suppose  that  by  some  mechanism  controlling  the  field 
particles,  the  number  of  field  particles  colliding  on  each  face  of  the  cube  is 
assigned  as  in  Case  1  - 1>^ »  and  further  each  particle  has  the  same  momentum 
components  as  given  by  (3.  18).  In  this  case,  <v>^  an<l  <^P  are  respec- 
tively  the  same  as  those  given  in  Case  1-b^.  However,  Ap  shows  no  fluctuation, 
and  <Ap  Ap>^  completely  vanishes. 
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2.  Multiple  Collisions  (t  j  0) 

The  conditions  are  similar  to  those  considered  in  Cases  1,  except  that 

the  test  body  and  the  field  particles  are  not  extremely  rigid  and  T  ,  the  period  of 

a  collision,  is  finite.  In  this  case  the  choice  of  At  is  limited  by  (2.  11)  and  <v> 

At 

given  by  (2.  13)  may  be  much  larger  than  unity: 

<v>At  =  £  wr  At  >  >  1  (2.  24) 

In  this  case,  Eq.  (2.  14)  is  still  valid 

<AP>At  ?r  Wr  At-  (2<14) 

This  is  the  same  as  what  was  obtained  in  Case  1.  For  <ApAp  >^  ,  however  , 
we  must  take  (2/26)  instead  of  (2. .22). 

< Ap Ap > t  =T—  I-*rr~  r*n  (p  +  —  +  p,_  )(p  +  ---  +Pp_)  (2.26) 

n  *  n 


Case  2-a .  First  we  consider  a  spherical  body  with  diameter  D  for  the 


test  body.  As  in  Case  1-a,  we  have 

1/2  ,  2  kT  ,1/2 


<V>(  .2* 


'n  ( 


m 


)  R2At, 


(3.20) 


_  4 

<Ap>  =  -  —  mv<v> 

*  At  3  At 

«-4  ^  —♦ 

For  calculating  <Ap  Ap>^  «  <Ap  Ap  Ap>^,  ....  it  is  convenient  to  use  Markoff's 
method  of  free  flights  as  briefly  explained  in  Appendix  D.  By  this  method,  as 
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shown  in  Appendix  E,  the  probability  of  Ap  being  between  P  and  P  +  dP  is  given  by 

1 


W(P)d  P  = 


x  exp 


_  3/2,  4  .  .  ,  _.3/2 

8  tt  (“<  v>At  mkT) 


+  0 


(  )?/2 
<V>At 


(3.21) 


“  <v  >  mkT 
3  At 


d3P 


See  Appendix  E,  Eq.  (E.  9).  By  means  of  W(P),  we  obtain 


<  Ap  Ap  >=  <Ap  Ap  >  =  <Ap  Ap  > 

x  xAt  y  y  At  z  z  At 


(3.22) 


*  f  mkT<V>At 


as  shown  in  Appendix  E.  So  far  as  <v>^,  ant*  <AP  ^P  are 

concerned,  there  is  no  difference  between  two  cases,  one  where  <v>^t  <  1 

and  the  other  <V>  >>  1.  However,  as  obtained  in  Appendix  E,  in  the  latter 

At 

<  Ap  Ap  Ap  Ap  >  is  proportional  to  <v>^t  .  It  is  not  simple  to  calculate  <Ap  Ap  Ap  >, 
since  this  involves  v,  and  hence  A(p)  is  not  uniform  with  respect  to  the  direction 
of  p  ,  In  the  following  cases, where  a  cube  is  taken  for  the  test  body,  we  calcu¬ 
late  those  quantities  more  easily. 

Case  2-b j .  In  order  to  make  treatments  easier,  we  consider  the  same 
test  body  as  in  Case  1  - b ;  that  is,  a  cube.  By  so  doing,  as  shown  in  Appendix  E, 

—4 

we  obtain,  by  assuming  v  =  (v,  0,  0), 
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and  again 


<Jp>t  =  -  j  mv<v>(. 


<Ap  Ap  >  =  <Ap  Ap  >  =  <Ap  Ap  >  =  f  mkT<v>  (3.23) 

x  x  At  y  y  At  z  rz  At  3  At, 


<Ap  Ap  >  =  <Ap  Ap  > 

x  y  At  y  z  At 


=  0 


<AP  Apx  Apx  Apx  >A?c  («V>  t)a 


(3.  24) 


So  far  as  <Ap>  .  and  <Ap  Ap>  ^  are  concerned,  there  is  no  difference  between 
At  At 

i. 

the  effect  of  binary  collision  and  the  effect  of  multiple  collision.  At  first 
glance,  one  might  be  puzzled.  We  note,  however,  that  in  these  two  preceding 
cases,  2-a,  2-bj,,the  distribution  of  field  particles  colliding  with  the  test  body 
is  at  random  twofold:  (i)  the  random  distribution  in  the  configuration  space  (on 
the  surface  of  the  test  body);  (ii)  the  random  distribution  in  the  momentum  space. 
It  is  possible  to  conceive  that  these  disorders,  or  one  of  them,  may  cause  such 
results.  In  order  to  investigate  the  situation,  we  calculate  <Ap  Ap>in  the 
following  three  cases. 

Case  2-b^.  On  each  of  the  faces  of  the  cube  which  are  parallel  to  v, 

<v>t 

— 7 -  field  particles  (with  no  fluctuation)  are  assigned  as  colliding  in  At;  on  the 

6  <v> 

A  t  V 

face  which  is  faced  in  the  direction  of  v,  ( — 7 —  +  —  D7  n  At)  field  particles*, 

o  Z 


♦These  numbers  are  the  same  as  in  Case  2-bj.  See  Appendix  E. 
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<  v> 

^  v 

are  assigned,  and  on  the  opposite  face,  (  — r—  -  —  DR  n  At)  field  partio  .  s  are 

o  2 

assigned.  The  distribution  in  the  configuration  space  is  then  no  longer 
at  random.  Assuming  that  those  particles  are  distributed  according  to  the 
Maxwell  function  in  the  momentum  space,  we  obtain  again 

<A?>t=  -  |  ®v<v>At 


< Ap  Ap  >  =  <Ap  Ap  >  =  <Ap  Ap  >  =-mkT<v> 

x  x  At  y  y  At  z  z  At  3  / 


At 


(3.25) 


<Ap  Ap  >  =  --  --  -  =  0 

x  y  At 

These  are  exactly  the  same  as  given  by  (3.  24). 

Case  2-bj.  We  now  assume  that  the  distribution  of  field  particles  in 
the  configuration  space  is  at  random.  However,  each  particle  is  assumed  to 
have  the  same  magnitudes  of  momentum  components 

IpxI  =  Ip  I  =  IpJ  =  JJJm  I  CJ  fd3c  '  JJJ  {d*c 


=  (mkT/2Ti) 


1/2 


(3.26) 


In  this  case  we  obtain  for  v  =  (v,  0,  0) 


<  v  >  = 

At  (2tt) 


6  nDn  .kT 

T72  <-> 


1/2 


At 


(3.  27) 


<AP>At  = 


'  3  m  V<V>At 


(3.  28) 
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<tp*  ‘"x  5t =  <£p,  ‘Vat  *  <4PZ  tp«  ^t 


=  A  mkT  <v> 

3  TT  At 


<Ap  Ap  >  = 

x  y  At 


=  0 


(3.29) 


Considering  the  relation  between  <Ap^t  and  <Ap  Ap  >  we  note  that  the 
temperature  of  the  test  body  in  equilibrium  is  no  longer  the  same  as  that  of 
the  gas.  This  conclusion  is  conceivable  since  we  ignored  the  Maxwell  distri¬ 
bution  of  particles  colliding  with  the  test  body. 

Case  2-b^.  Finally  we  assume  that  the  number  of  particles  colliding 
on  each  face  of  the  cube  is  assigned  as  in  Case  2-b^,  and  further  each  particle 
has  the  same  momentum  component  as  giyen  by  (3.  26).  In  this  case  <V^ 
and  <AP>At  are  respectively  the  same  as  those  given  in  Case  2-b^.  However, 
Aphas  no  fluctuation,  and  <Ap  Ap>t  completely  vanishes. 

These  results,  obtained  so  far  through  rather  primitive  (but  precise) 
treatments,  reveal  some  aspects  of  multiple  interaction.  So  far  as  <v>^t  > 

-4 

<Ap>  and  <Ap  Ap  >  are  concerned,  there  is  no  difference  in  effect  between 
At  At 

binary  collision  and  multiple  collision.  The  binary  interaction  is  based  on  the 
hypothesis  of  molecular  disorder.  Our  treatment  of  multiple  interaction  is  also 
carried  out  by  granting  the  same  assumption.  The  assumption  is  justified  by 
assuming  that  mutual  interactions  among  field  particles  are  negligible.  As  a 
matter  of  physics,  the  binary  collision  assumption  and  the  assumption  of  no 
mutual  interaction  among  field  particles  are  compatible.  It  is  noted,  however, 
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that  the  multiple  collision  assumption  and  the  assumption  of  no  mutual  inter¬ 
action  among  field  particles  are  not  necessarily  compatible,  as  considered 
in  the  following  section. 

Another  important  aspect  is  that  <Ap>At>  (and  also  <ApApAp>At),  is 
directly  related  to  the  non-uniform  distribution  of  field  particles  reflected 
from  the  surface  of  the  test  body:  In  other  words,  <Ap>At  is  said  to  be  caused 
by  the  wake  produced  behind  the  test  body.  On  the  other  hand,  the  wake  has 
no  effect  on  <ApAp>At,  so  far  as  mv  /kT  <  <1:  <ApAp>At  is  caused  simply  by 
fluctuations  which  appear  in  Ap,  and  are  smoothed  out  in  <Ap>A^«  See  the 
details  of  calculation  in  Appendix  E: 

SECTION  IV 

AN  ELECTRICALLY  CHARGED  TEST  BODY  IN  ELECTRON  GAS 

In  Section  III,  we  obtained  experimentally  a  rule  that  multiple  inter¬ 
actions  are  equivalent  to  superposed  binary  interactions,  when  mutual  inter¬ 
actions  among  field  particles  are  ignored  and  so  far  as  the  friction  and  the 
diffusion  (in  the  momentum  space)  of  the  test  body  are  concerned.  In  this 
section,  we  assume  that  the  rule  is  valid  for  considering  interactions  among 
charged  particles. 

We  consider  for  the  test  body  a  particle  with  electric  charge  Q  and  mass 
M.  The  test  body  is  submerged  in  an  electron  gas;  a  positive  charge  is 
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assumed  to  be  continuously  and  uniformly  spread  in  the  space  so  that  the 
positive  charge  neutralizes  the  charge  of  the  whole  electrons.  We  also 
assume 

M  >>  electronic  mass  m  (i) 

Q  <  electronic  charge  e  (ii) 

and 

mv2/kT  <<1  (iii) 

where  T  is  the  temperature  of  the  electron  gas  which  is  in  thermal  equilibrium, 
v  the  speed  of  the  test  body  relative  to  the  electron  gas. 

As  stated  in  the  Introduction,  the  interaction  among  those  particles  is 
due  to  (a)  the  forces  exerted  on  the  test  body  by  field  particles  (electrons),  (b) 
the  forces  exerted  on  field  particles  by  the  test  body,  (c)  forces  among  field  particle* 
Approximation  A 

Let  us  first  ignore  (b)  and  (c)  and  see  the  consequence  of  applying  the 

rule  set  forth  at  the  beginning  of  this  section. 

Effect  1.  <Ap  >  =0.  Consider  the  effect  of  an  electron,  with 

-  r  At 

velocity  c  which  is  parallel  to  v,  encountering  the  test  body  with  impact  vector 
r.  [Here  impact  vector  is  defined  as  impact  parameter  together  with  its 
direction,  ]  We  may  consider  another  electron  with  the  same  velocity  with  impact 
vector  -r*.  The  effects  cancel  each  other.  Next  we  consider  an  electron,  with 
velocity  c  which  is  perpendicular  to  v,  passing  by  the  test  body  with  impact 
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parameter  r.  We  may  consider  another  electron  with  the  same  velocity  and  with 
impact  parameter  -r.  The  effects  of  theae  two  electrons  cancel  each  other. 

See  Fig.  2.  Considering 


Fig.  2.  vj  the  velocity  of  the  test  body 
c:  the  velocity  of  two  electrons 
r:  the  impact  vector. 

The  momentum  given  to  the  test  particle  by  the  two  electrons 
vanishes.  Here  the  electrons  are  free  from  any  of  the  other 
particles. 


the  density  of  electrons  dependent  only  on  Jc  |,  we  may  conclude 

<  Ap>^t  =  0  (4. 1) 

Effect  2.  <Ap  Ap>^t  =  oo.  An  electron  with  velocity  c-v  relative  to  the 

test  body  and  impact  parameter  r  gives  the  test  body  momentum 

-  f+°°  eQr  dt 

Pi  =  J  -oo  "tri  +  f/2 

=  -r le-9—  l 

r  |  c- v  |  r 


A  test  body  suffers  a  deceleration  drag  in  average  when  it  moves  through  force 
fields  distributed  at  random,  because  the  total  period  of  deceleration  is  longer  than 
the  total  period  of  acceleration.  The  average  drag  is  of  the  order  of  F2£/ (Mv2), 
where  F  is  the  average  intensity  of  the  force  fields,  1  the  correlation  length  of  the 
fields,  M  the  mass  of  the  test  body,  v  the  velocity  of  the  body.  It  is  assumed  that 
FJf<  <Mv®.  This  sort  of  drag  is  considered  small  and  is  ignored. 
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Hence 


_  _  4epCy 
P  iPj  =  rac2 


(4.2) 


2,  2 


by  ignoring  v  /c  according  to  (iii).  Remembering  Cases  1-b  and  2-b  in  the  last 
section*,  and  considering  interactions  as  if  binary,  we  may  calculate  <ApAp>t> 
by  taking  f  for  the  distribution  function  of  electrons. 

<A?Ap>t  =  JJ]J  7T“F-  2n  rdrfcd3cAt 


8ireaCy*  At  JjJ  fdacJ^ 


(4.  3) 


Here 


m  .3/2 


—  d3c  =  n  (  ) 

c  2nkT 


JIT  <- 


m 

2kT 


c8)  sin  0  d  9  dcp 


/  m  .3/ 2  „  2kT 
n  2nkT  4"~ 


4n  .  m  1/2 
JW  {  2kT  ' 


(4.  4) 


J  ±'  -  [log  .f- 


*In  the  Approximation  of  ignoring  mv2/kT,  we  do  not  need  to  consider  v  in  the 
calculation  of  <Ap  Ap>^t  . 
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The  result  diverges  either  for  rj  =  0  or  for  Xg  =  oo. 

Approximation  B.  Secondly,  we  ignore  (c),  mutual  interaction  among 

— *  “4 

electrons,  and  take  into  account  (a)  and  (b).  Here  <  A  p  >.  .  and  <ip4p>,.  are 

Qt  At 

calculated  by  taking  all  the  interactions  as  binary  and  the  results  are  well-known: 

lim  <  A  p  >  =  0  [  log  rs  ]  =  oo, 

rs  -*  oo 

(4.  5> 

lim  <  A  p  A  p  >  =  0  [  log  ra  ]  =  oo 

In  the  two  cases  stated  above,  spatial  fluctuations  in  the  electron 
distribution  are  permitted  with  no  restriction:  Even  the  probability  of  all  the 
field  electrons  to  come  together  at  one  spot  is  taken  into  account.  On  the  other 
hand,  the  interaction  of  the  Vlasov  type  is  conceived  by  assuming  that  there  is 


*  Fluctuations  ol  the  electric  force  field  in  an  ionized  gas  were  calculated  by  the 
Markoff  methof  of  random  flights  by  J.  Holtzmark  [Ann.  d.  Physik  58,  577(1919): 
Physik,  Zeits.  20,  162(1919),  25,  73  (1924)  ].  Later,  Chandrasekhar  and  von 
Neumann  \  4  ..  used  the  same  method  for  calculating  fluctuating  forces  exerted  on  a 
star  by  other  stars.  They  considered  not  only  the  force  but  also  its  time  derivative 
due  to  an  assumed  distribution  of  star  velocities.  In  the  initial  formulation,  they 
considered  the  correlation  of  the  two  quantities  (force  and  its  time  derivative)  of 
each  field  star.  Because  of  mathematical  difficulties,  they  abandoned  the  precise 
correlation  between  them  and  made  the  average  of  one  at  a  given  value  of  the  other. 
Then  they  calculated  the  correlation  period  of  a  force  by 

(a  given  value  of  force)/(the  average  time  derivative  of  force  at  the 
relevant  value  of  force.  ) 

The  result  apparently  converges.  We  notice,  however,  that  the  convergence  of  the 
result  is  not  proved  by  such  an  approximate  treatment  as  done  by  Chandrasekhar 
et  al. 
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no  spatial  fluctuation  in  the  electron  distribution.  The  real  (feasible)  condition 
must  be  between  those  two  extreme  conditions.  In  the  following,  we  consider  the 
condition.  This  is  somehow  similar  to  a  condition  of  turbulence,  in  which  dynamic 
characteristics  of  a  fluid  set  restrictions  to  fluctuations  appearing  in  the  flow  of 
the  fluid.  ^ 

Approximation  C.  Finally,  we  conclude  that  the  divergences  (4.5)  are  attributed 
to  the  neglect  of  mutual  interactions  among  electrons.  The  widely  accepted 
solution  is  to  consider  the  polarization  in  the  distribution  of  electrons  due  to  the 
potential  field  induced  by  each  of  the  particles  including  the  test  body.  As  stated 
in  Part  I,  in  general,  this  solution  does  not  seem  plausible.  Then,  what  is  the 
effect  of  mutual  interactions  among  field  electrons  by  which  the  result  converges? 
We  note  that,  by  the  Baltzmann  type  binary  interaction,  regardless  of  the  type 
of  potential  between  two  interacting  particles,  a  limitless  fluctuation  in  the  spatial 
distribution  of  field  particles  is  permitted  to  appear.  In  other  words,  each  of  the 
field  particles  interacting  with  the  test  particle  is  permitted  to  appear  anywhere, 
regardless  of  the  distribution  of  the  other  field  particles.  Our  assertion  is  that 
such  limitless  fluctuations  in  the  spatial  distribution  of  field  particles  in  inter¬ 
action  with  the  test  particles  are  not  permissible  in  an  electron  gas,  since 
potential  energies  among  field  electrons  prevent  those  electrons  being  accumulated 
at  local  spots.  Such  fluctuations  may  be  permitted  to  appear  only  in  the  distribution 
of  neighboring  electrons  of  the  test  particle. 
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Fig.  3.  0:  the  test  body,  1,  2,  3, - :  trajectories  of  electrons. 

If  the  velocity  of  the  test  body  with  a  large  mass  is  zero, 
each  electron  which  passes  by  near  the  test  body  is  rejected 
toward  the  outside,  and  anomalies  occur  in  the  distribution 
of  electrons.  Such  a  anomalous  distribution  near  the  test 
body  is  not  Maxwellian. 


Considering  these,  we  propose  the  following  scheme  of  treating  electrons 
in  an  electron  gas: 

1.  The  interaction  of  a  test  electron  with  its  mutual  nearest  neighbor 
is  treated  as  binary  interaction. 

2.  When  observed  by  an  observer  resting  on  the  test  electron,  the 
distribution  of  the  mutual  nearest  neighbors,  which  come  to 
interaction  with  the  test  electron  one  by  one,  is  not  uniform  in 
general:  If  the  test  electron  has  no  velocity  relative  to  the  average 
velocity  of  the  other  electrons,  the  density  distribution  of  the  nearest 
neighbors  is  lower  with  a  spherical  symmetry  surrounding  the  test 
electron.  If  the  test  electron  has  a  velocity  relative  to  the  average 
velocity  of  the  other  electrons,  the  density  distribution  to  the  nearest 
neighbors  is  lower  behind  the  motion  of  the  test  electron.  In  other 
words,  there  is  a  wake  behind  the  motion  of  the  test  electron. 
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3.  Each  of  the  electrons  may  assume  a  test  electron.  Thus  each  electron 
has  its  own  wake. 

4.  Except  for  the  wakes,  all  the  electrons  are  seemingly  uniformly  distribute 
with  no  fluctuation  in  a  similar  manner  as  considered  by  Vlasov. 

5.  An  electron  accompanied  by  its  mutual  nearest  neighbor  exerts  a  fluctuat¬ 
ing  force  on  other  electrons  because  of  the  anomolous  distribution  of  its 
nearest  neighbor. 

The  kinetic  theoretical  treatment  of  an  electron  gas  according  to  this  model 
will  be  formulated  rather  schematically  in  the  following  section. 


SECTION  V 

A  SCHEME  OF  KINETIC  THEORETICAL  TREATMENT 

* 

OF  AN  ELECTRON  GAS 


It  is  easily  shown  that  the  Liouville  equation  is  reduced  to 


dF^x.  ;t)  p.  dFxfx.jt) 


at 


m 


+  S  Vs  y -~l 

j-*--k 


*xk:t) 


dx.-  -  -dx,  =  0 
J  k 


(5.  1) 


This  schematic  presentation  lacks  detailed  correlation  functions  and  is  not 
suitable  for  precise  analyses.  See  Part  III. 
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i  #  j  * - ♦  k. 


We  define  W  as  the  probability  of  two  particles,  i  and  j,  being 
mutual  nearest  neighbors:  W.^  is  a  function  of  (q.  -  q  )  and  of  the  number 
density  of  electrons. 

By  putting 


F2  (ij)  *Fa  (x.x^;t), 


F(ij,  =  Wij  F2  (ij>' 


(5.2) 


F i/ j  »  O-Wy)  F*  (ij) 


We  have 


F,  (ij)  =  F(ij)  +  F.,. 


(5.3) 


Since  the  interaction  between  i  and  j  is  considered  weak  in  the  domain  where 

Fi/j  Predominate8-  we  maY  Put  Fi/j  =  tFi(i)  Fi  (j)  +  Fa(ij)  ] 

where 

FjU)  2FX  (x.;t) 

By  ignoring  Fa  (ij),  we  have 

Fi/j  =  (1"Wij)Fi(i)  Fi(j)  (5.4) 
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Hence  we  write  for  (5.  3) 

Fa(ij)  =  F i(i)  F i(j)  +  Wtj  [Fa(ij)  -Fx(i)  Fx(j)] 


With  respect  to  Fa(x.x.x,  ;t),  we  similarly  define 

1  J  K 

Fs  (ijk)  =F3  (x.xjxk;t), 


F(i,)/k  *  w.j  F»  <«“>• 
F(ik)/j  ■  wik  F=  <«“>■ 
F(jk)/l  ■  wjk  F»  ™k> 


F...  ..  =  (l-w..  -W..  -W.,  )Fs  (ijk) 
i/j/k  lj  lk  jk'  w 


In  the  same  approximation  as  of  (5.  4),  we  write 
F(ij)/k=  Wij  Fs  (ij)  F>  (k)’ 

F(ik)/j  "  Wik  F’(ik>  F>  «>’ 


F0k)/i  =  wjk  F»«k)  F>{i) 


Fi/j/k  =  (1-  wij-wjk-  wik)F1(i)  F>(j)F1(k)- 


F3  (ijk)  =  Fx(i)  Fx(j)  F x(k) 

+  W..  [  Fa(ij)  -  Fj(i)  Fx(j)  ]  Fx  (k) 
+  W.k  C  F s(ik)  -F |(i)  Fi(k)  ]  Fx  (j) 


(5.5) 


(5.6) 


(5.7) 


(5.8) 


=Fl(i)  Fj(j)  F x(k) 

+F 2  <ij)  Fx(k)  +  Fa  (ik)  FiO) 
+F a  (jk)  Fx(i) 


(5.9) 


obtain 


where 


If  we  take  3  for  s  in  Eq.  (5.  1)  and  substitute  (5.  9)  in  the  equation,  we 


SF^i)  p  JF,(i) 

+  — +  i|+fe+$3+^*=0 
m  oq. 

l 


at 


(5.  10) 


*1  = 


±  y 

V  /. 


J  «- 


J  FiU)  dx. 


•  F>  (i) 
api 


(the  effect  of  Vlasov  type 
force) 


v-  1 1  JI { ^ij  +  yik>  F*  (jk)  dxj  dxk  ’  W-  Fl  (i) 

j  k 


(5.  11) 


(5.  11). 


(the  effect  of  fluctuations  in  the 
field  particle  distributions) 


=■  V  i  J  ?  •  r'  «>  dxj 


IJ 


(5.  ID, 


1  1 


(the  effect  of  Boltzmann  type 
interaction  among  mutual 
nearest  neighbors) 

a 


** =  t  v*  L  L  J  J  -7ik  Fl  <k)  dxk  •  af.  F*  (ij)  dxj 


j  k 


(5.HK 


(the  effect  of  Vlasov  type  force 
on  the  correlation  between 
mutual  nearest  neighbors). 
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By  taking  4,  5, - for  s  in  (5. 1)  we  may  have  more  complex  effects.  The 

detail  will  be  investigated  in  Part  III  of  this  report, 

SECTION  VI 

CONCLUDING  REMARKS 

1.  According  to  the  investigations  carried  out  in  Sections  II  and  III, 

it  is  most  likely  that  the  effect  of  multiple  interactions  of  a  test  body  with  field 
particles  is  the  same  as  the  effect  of  interactions  as  assumed  to  be  binary, 
so  far  as  the  friction  and  diffusion  of  the  test  body  in  momentum  space  is 
concerned,  and  so  far  as  the  mutual  interactions  among  field  particles  are 
ignored.  The  difference  between  two  modes  of  interaction,  multiple  and 
binary,  appears  when  the  mutual  interactions  among  field  particles  are  taken 
into  account.  We  induce  this  conclusion  from  the  results  of  our  various  mental 
experiments.  We  might  not  be  allowed  to  claim  the  conclusion  to  be  a  general 
law.  At  least  we  may  be  allowed  to  propose  the  conclusion  as  a  hypothesis. 

2.  In  view  of  the  above  conclusion,  we  consider  an  electron  gas.  The 
mutual  interactions  among  field  electrons  prevent  them  from  being  in  free 
flight.  Effective  fluctuations  in  the  distribution  of  field  electrons  are:  i) 
the  fluctuations  in  the  distribution  of  mutual  nearest  neighbors  of  the  test 
electron  caused  by  the  test  electron;  ii)  the  fluctuation  of  each  field  electron 
caused  by  its  mutual  nearest  neighbors. 

3.  A  scheme  of  kinetic  theoretical  treatment  of  an  electron  gas  based 
on  the  Lionville  equation  is  proposed.  There,  the  assertion  (2)  stated  above 
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is  formulated. 


SECTION  VII 
REFERENCES 

1.  Einstein,  A.:  Investigations  on  the  Theory  of  the  Brownian  Movements. 

Edited  by  R.  Fiirth,  Dover  Publications,  Inc.  ,  1956. 

2.  Uhlenbeck,  G.  E.  and  Ornstein,  L.  S.  :  Phys.  Rev.  36,  823  (1930). 

3.  Wang.M.  C.  and  Uhlenbeck,  G.  E. :  Rev.  Modern  Phys.  ,  17,  323  (1945). 

4.  Chandrasekhar,  S. :  Rev.  Modern  Phys.,  J_5,  1  (1943). 

5.  Wax,  N.  (Editor):  Selected  Papers  on  Noise  and  Stochastic  Processes  by 

Doob,  Ornstein,  Uhlenbeck,  Rice,  Kac,  Wang  and  Chandrasekhar. 
Dover  Publications,  Inc.,  New  York,  1954. 

6.  Batchelor,  G.  K.  :  The  Theory  of  Homogeneous  Turbulence.  Cambridge 

University  Press,  1956. 

7.  Bierens  de  Haan,  D.  :  Nouvelles  Tables  D'Integrales  Definites.  Hafner 

Publishing  Company,  New  York,  1957,  pp.  394,  506. 

8.  Jeans,  J.  H.-:  The  Dynamical  Theory  of  Gpses.  Dover  Publications,  Inc.  , 

New  York,  Fourth  Edition,  Appendix  A. 


36 


APPENDIX  A 


THE  DERIVATION  OF  THE  FOKKER- PLANCK  EQUATION 
FROM  THE  SMOLUCHOWSKI  EQUATION 

By  assuming  |  Ap|  <  |p|  and  expanding  functions  in  (2.  1)  in  Taylor's 
series,  we  have 

_  af  ,  aaf  aaf 

yp.t+At)  =  fo(p,t)  +  -gj-At  +  j-,  -^(Atp+3-|  +.....  , 


_  Sf  .  asf 

fo(p  -  At,  t)  =  fQ  (p,t)  -  Ap  .  +Ti  Ap  Ap  :  ^ 

..  -  a»f 


3!  Ap  Ap  :  +  *’ 


1)1  (p  ’.t  |p,  t  +  At)  dap' 


•'[«(;,  t;  ip,  it) -if  .  H-+2Ts?4?- 

Jirfe +  ~  ]d’4p 

as  given  by  (2.  4).  On  substituting  these  in  Eq.  (2.  1),  we  obtain 


f  (P,t)  +  at  +jl  i!f° 

°  at  zj  at 


I  2  (  At)  S  + 


=  JJJ  C  (fo  •  Ap 


at 


o  ,  1 


a2f 


a?  ’'***■  m ? 


^3  f 

-  JT  Ap  Ap  Ap  :  >--)  ®(p.t;Ap,  At) 
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+  (fO 


-  Ap 


1  .  -  .  - 


+  21 Ap  Ap:  W? 


-) 


a£a£  . 
2  ! 


3  *co 

dp  dp 


....  ] d3Ap 


(A-l) 


According  to  (2.  5)  and  considering  that  Ap  is  independent  of  p,  we  have 

Hj  ^  vi‘t‘p 


5?  •  <4p=a, 


(A-2) 


<4pa>’t 


Hence  (A-l)  yields 


df  33f 

*T  At+2-  17T(^,3 


df 


d3f 


d3f 


<Ap>At  •  Tf  +  fi  <ApAp>At:dT5?  '?T<ApApAp>At:dFd?d? 


7T  !%k  k  :<ApApAp>At  + 


+  2  ! 


3s 


f  3  f 


+  f  -S— :  :  <Ap  Ap  >  -  ^ 

odpdp  ^  ^ 


o\  * 


4t"y»»  Jlpl f  :  <'p''p*p’m  + 
-T.[—  J+— - 


^  •  (</iP>&tfo,+  2!T^T:  (<^P^P>Atfo) 


-  -si "ITap a I?  :  {  A?  A?  >At  V 


+  I 


— 2 -  :  (<&pAp&p&p>  f  )- 

4  !  dp  dp  dp  ap  •  At  o 


(A.  3) 


If  it  is  shown  that 


<  Ap  >^/At  =  <AP  > 


<Ap  Ap  5^  /At  =  <  A  p  Ap  > 


are  independent  of  At  and  that 


dt*  ^  ^ ^  At  <  < 


0(4.)=/(f£«)<  <1 


(A.  4) 


(A-4) 


under  certain  conditions  imposed  by  physical  conditions  of  our  interest, 
Eq.  (A-3)  is  the  Fokker-Planck  equation. 


39 


APPENDIX  B 


THE  FOKKER- PLANCK  EQUATION 
DERIVED  FROM  THE  BOLTZMANN  EQUATION 

The  Boltzmann  equation  is  based  on  the  assumption  of  binary  collision. 
The  collisions  are  Markoff  processes.  Hence,  if  most  collisions  are  weak,  both 
regarding  the  test  body  and  the  field  particles,  the  Boltzmann  equation  must  be 
reduced  directly  to  the  Fokker-Planck  equation  in  the  form  as  given  in  Appendix  A. 
Usually,  *  the  derivation  is  carried  out  by  means  of  equations  of  moments  based 
on  the  Boltzmann  equation.  Here  is  given  a  direct  derivation:  The  Boltzmann 
equation  by  means  of  the  conventional  notations  is  written  as  follows: 

df 

dT=  III  Cfl  >  fo  (P'>  ‘  f>  (Pj)  lo  ^  ] 

x  B  db  d  e  d^j 

Here  ft  is  the  distribution  function  of  the  field  particles,  b  the  impact  para¬ 
meter  of  collision  between  the  test  body  and  field  particles,  e  the  longitudinal 
angle  of  the  plane  on  which  the  trajectory  of  a  colliding  molecule  is  present, 
and  B  a  function  of  b,  the  relative  velocity  and  the  force  between  the  test  body 

and  the  colliding  molecule.  There  are  certain  relations  among  p/,  p' ,  px  and 

— •  _ 

p:  Suppose  that  the  collision  between  the  test  body  with  p  and  a  field  molecule 
with  px  under  the  condition  of  impact  parameter  b  causes  the  test  body  with  p' 

*W.P.  Allis:  Motion  of  Ions  and  Electrons,  in  Handbuch  der  Physik,  edited  by 
S.  FKigge  (Springer,  Berlin,  1956),  Vol.  21,  P.430 
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-*4  f  ^  f 

and  the  field  molecule  with  px .  Then  the  collision  between  the  test  body  with  p 

—4  f  • 

and  the  field  molecule  with  plunder  the  same  condition  of  impact  parameter  results  4a 

■4  -4 

the  test  body  with  p  and  the  field  molecule  with  px.  This  relation  is  shown  by 
considering  that  the  collision  process,  according  to  Newton's  equations,  is 


Fig.  4.  In  (a),  a  particle  with  momentum  p  and  another  particle  with 
momentum  px  collide  and  their  momenta  after  the  collision 
are  p '  and  px  .  The  relation  is  reversible  as  shown  in(b). 

By  putting 

Ap  =  px'  -  pi  =  p  -p'  (B-2) 

Ap  is  a  function  of  Pi.  P.  b.  and  hence  we  may  take  for  the  independent 
variables 

Ap,  p,  b,  0  (B-3) 

instead  of 

Pi.  P,  b,  fl; 
or  we  may  write 

Pi  =  Pi  (  Ap,  p,  b,  9).  (B-4) 
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By  considering  the  reversible  relations  between  (plt  p)  and  (p^.  p*),  we  may 
say  the  following:  If  we  define  f*i  by 


f*i  (p;  Ap)  =  ft  (pi), 

then,  considering  (B-2),  (B-4),  we  have 

(pi)  =  fi*  (p -  Ap)  =  f1*(p-Ap:  -Ap) 
By  expansion,  we  obtain 

Mpi  )  =  fi*(p;-Ap)  -  Ap.  ^ 

d  p 


(B-5) 


(B-6) 


,  i  a’f,*  1_  ,  v  , 

+  27ApAp:  3^  *  31  (  )+" 


(B-7) 


at 


a3f 


fo(p'Ap)  =  £o(p)  • &p  •  aT  +  z7  Ap  Ap  :  + 


(B-8) 


Further  we  define  B*  by 

d  (pyc.Px  y,  Piz} 
B*  =  B  .  9  (Tpx,  Apy,  Apz 


(B-9) 


Equation  (B-l)  yields 
d  f 

IT  =  ID [f?(p- Ap:  Ap  )  {o  (P  '  Ap^ 

(B- 10) 

-fj*  (p;  Ap)  fQ  (p  )  ]  B*  bdbde  d3  Ap 
Substituting  (B-7),  (B-8),  and  (B-9)  in  (B-10),  we  obtain 
d  f 
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df 


Sfi*Mp  )  +  {jfl  (_Ap )  Ap.  ] 


a? 


df 


♦  ^  Cfc  *  * :  i^£l  +  2  A?  A> !  ^ 


aaf 


+  fi*(-  Ap  )  Ap  Ap  : 


]-  r.  C  ] 


apap  3! 


+ — 


B*  bdbded3  Ap 


(B-ll) 


Considering  that 

,.+oo 


-oo 


ff*  (-  Ap)  d3  Ap=  JJJ  +<jo-f t  (  Ap)d3Ap 
=  JjJ+°°  fi*  ( A  p)d3  Ap 


we  have 


JJJ  [fj*(-Ap)  fQ  -  fx*(  +  Ap)  fQ  ]  B*  bdbded3 A p  =  0 

Further  we  define 

JJJ  Ap  f**  (-Ap)  B*  bdbded3Ap  =  <  Ap> , 
jJJ  Ap  Ap  f^(-Ap)  B*  bdbded  3Ap  =  <Ap  Ap> 


(B- 12) 


Equation  (B-ll)  now  yields 


df 
_ o 

dt 


a 

2p 


<Ap>fo  + 


_1_ 

2! 


— —  :  <Ap  Ap  >f  - 
apap  o 


]+  (B-13) 
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APPENDIX  C 


A  SOLID  AND  ELASTIC  TEST  BODY  IN  A  RAREFIED  GAS 
CONSTITUTED  OF  SOLID  AND  ELASTIC  MOLECULES 


I.  A  Spherical  Test  Body  (T  ig  almost  zero) 

We  calculate  <Ap>,  <Ap  Ap>,  . . .  by  assuming  that  the  test  body  and 
the  molecules  constituting  the  gas  are  extremely  hard.  Due  to  the  assumption, 
each  collision  period  is  extremely  short  and  we  may  take  At  as  short  as  we  like. 
Consequently  the  collisions  are  binary. 

First  the  test  body  is  assumed  to  be  a  sphere  of  radius  D/2,  mass  M 
and  velocity  relative  to  the  gas  v,  while  the  molecules  are  spheres  of  radius 
a/2  and  mass  m.  Further  the  gas  is  assumed  to  be  in  thermal  equilibrium  with 
temperature  T.  Here 


M  >  >m 


kT  * 

v  <<  (^) 
m 


(C-l) 

(C-2) 


are  assumed.  By  taking  a  spherical  coordinate  system  with  origin  O  fixed  at 
the  center  of  the  test  body,  the  coordinates  of  the  center  of  a  molecule  colliding 
on  the  surface  of  the  test  body  are  (R,  9,  cc),  where 


R  =  (D  +  o)/2 

9  the  colatitude,  and  cp  the  longitude.  The  number  of  molecules  which  have 


velocities  between  c  and  c  +  dc  is  given  by 


fd3c  =  n  ( 


m 


2  nkT 


3/2  r  m 
)  exp  L 


2+cH+cs)] 

2  k  T  x  y  z 


d3c 


(C-3) 
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where  c  is  measured  by  an  observer  at  rest  with  respect  to  the  laboratory 
(not  to  the  test  body  which  is  moving  with  velocity  v).  The  number  of  such 
molecules  which  collide  with  the  test  body  on  an  elementary  surface  area 


dS  =  Ra  sin  0  d  0  dcp 

during  the  period  between  t  and  t  +  At  is  given  by 


(C-4) 


where 


A  n  =  fd3c  |  V|  cos  0  dS  At 
cs 


(C-5) 


— ♦ 

V  = 


c-v 


(C-6) 


and  the  z-axis  (polar  axis)  is  chosen  so  that  the  axis  is  parallel  to  V  in 
the  opposite  direction. 


Fig.  5.  ds  =  R*  sin  0  d0dcp 

By  carrying  out  the  integration  with  respect  to  0,  from  0  to  tt/2,  and  cpfrom 
0  to  2tt,  we  obtain 


A  n  =  ttR®  I  vl  fd3c  At  (C-7) 

c 

By  choosing  a  rectangular  coordinate  system  (?, r|,  C  )  so  that  5-axis  is  in  the 

—4 

direction  of  v,  we  have 


45 


V  =  (v,  o,  o). 


c.  -  v  =  C_ , 

?  s 


c  =  C  , 
n  T) 


V  Cc 


{ ■  ”  < SIt)3/2  -  ITt  [  (cs  +  v,a  +  V  +  ci  : 


(C-8) 


dac  =  d3C 


V?  =  C.a  +  C«  +  Cf- 

?  n  C 


Substituting  these  in  (C-7),  we  may  carry  out  the  integration  with  respect  to  c  . 


A  n  =  ttR3  At  JJJ  Cfd3C 


(C-9) 


By  considering  (C-2),  we  have 


exp[-  ^(C§+v)3] 


=  exp  ['^rcp2  ]x  [1  •  S(2c§v  +  vS) 


+  rr  <2kTr  (2C? v+v2)2  *  r. (  )  +  --] 


(C- 10) 


=  eXP[-^  C/  ]d-fT2C5v) 


Hence 


6n--nR3„j4n„(  f  C*  «q>  (-  C‘,  dC 


+  00 


-f  ft  CC*  "■>  [  -  2Tt  (CfS*  Cn  +  cc’  >  1  d3c 


-oo 


,  1/2  .  2k T  ,1/2  A  + 

2  n  n  R-  (  -  )  At 

m 


(C-ll) 
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(i)  <Ap>  .  Each  of  those  molecules  gives  the  test  body  momentum 


6p  in  the  direction  of  R: 


6p  =  -  2  m  j  V  |  cos  0  S./R 


(C-12) 


(6p  )  =  6  p  sin  0  sin  cp, 
(6p)  =  6p  sin  0  cos  cp, 

y 

(6p)  =  6  p  cos  0 . 

z 


(C- 13) 


By  giving  the  direction  of  v  in  the  present  coordinate  system  by  0  =  0^,  cp=  0, 


we  have 


cos  0  = 
o 


v3  +  V3  -  cs 


V3  -  V.  c 

IvTlvl 


(C- 14) 


and  the  component  of  6p  in  the  direction  of  v  is 


(6p)  =  (6p)  sin  0  +  (6p)  cos  0 

xiv  y  o  z  o 


2m  | V |  (cos3  0cos  0q  +  sin  0  cos  0  cob  co  sin  0Q  )  (C- 15) 


The  momentum  transferred  to  the  test  body  from  the  field  molecules  with 

velocities  between  c  and  c  +  dc  in  At  has  its  component  in  the  v* direction 
2n  „Tt/2 

A  p  =  i  -  2m  |v|  (cos3  0  cos  0  +  sin  0  cos  0  cos  cp  sin  0  ) 

c  liv  J  o  o 

'  cp=o  0  =  0 


.  x  f  | V |  R*  cos  0  sin  0  dfidcp  d3c  At 


=  -  tt  mRs  Vs  cos  0  f  d3  c  At 
o 


=  -  n  mRa 


-v.  c)  lc-v 


fd3c  At 


(C- 16) 
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For  carrying  out  the  integration  with  respect  to  c,  we  remember  (C-8)  and 
obtain 


^-TTmR-AtjJj  -  CFC(l-  c?v) 


i  m  %3/2 

x  n  (  mt!  exp 


m 


2kT  (c?2  +  cn2  +  cC  »  *  d3° 


=  -  tt mR2  At  n  (  )3/2  U1 

2nkT  kT 


X  IfJ  C5»Cve!ip(--i|-TC-,d=C 


By  putting 


C 


§ 


C  sin  9X  sin  CO i 


daC  =  C  r’dC  sin  0!  d  81  dcpi 


We  ha  ve 

jffc.SC«xp(-  ^C-ld'C 
=  [JJ  sin3  9j  sinncp1  d9x  dcpj  ] 

x  [  j  C-  exp  (  -  ^  Cs)  dC] 

-  4  n  f  2kT 

3  '  m  ’ 

Hence 

A  Hi  v  =  -  -|-  n1^2  mR'  n  (  v  At  +  0  [  v2] 


(C-17) 


The  component  perpendicular  to  v  does  obviously  vanish  because  of  the  axial 
symmetry  of  the  conditions  with  respect  to  the  direction  of  v.  According  to 
to  the  notation  given  in  Section  II,  we  may  write 
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(C- 18) 


<&^t  = 


-  —  m  v  An 
3 


(ii)<ApAp^.  We  first  calculate 


<{Ap)a>=  <Apx  Apx>+  <Apy  Apy>+  <Apz  Apz>  (C-19) 


According  to  (C-5)  we  have  for  this  part  of  <Ap  Ap>  due  to  the  molecules  of 
velocities  between  c  and  c  +  dc 

<  (Ap)8  >  =  At  fd3c  J  (6p)2  | V |  cos  0  d  S 
s 

where  -6pis  given  by  (C-12).  By  putting 


We  obtain 


dS=R8  sin  0  d^dcp 


<  (Ap)8>  =  2  rtm2R8  V3  fd3C  At 


or  by  considering  (C-10) 


<(Ap)5,>At  =  8  tt8R-  n  m2  (~)  3/2At  +  0  [v8] . 


m 


(C-20) 


Next,  we  calculate  <Ap1iV  8  >with  respect  to  the  component  Apuv  in  the  v  - 

direction.  By  considering  (C-15)  and  (C-10),  it  is  easily  shown  that 

<(Apilv  )8  >  =  ■j  <  (  Ap)8  >  +  0  [  v?]  (C-21) 

Therefore,  we  may  conclude  that  <Ap  Ap  >.  <Ap  Ap  >,  and  <Ap  Ap  >  are  inde- 

x  rx  At  y  ry  At  z  *z  At 

pendent  of  v ,  and  4 

<Ap  Ap  >  =  <Ap  Ap  >A  =  <Ap  Ap  >,  =  T  m3  An  (  — — ) 

x  At  y  *y  At  *z  rz  At  3  nn 
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The  other  components  are  of  the  order  of  v®,  and  hence  are  ignored. 

•4^-4  ^ 

By  means  of  similar  ways,  we  may  obtain  <Ap  Ap  Ap^t>  ^pApApAp^, .  . . 
An  essential  and  common  feature  of  those  results  is  that  they  are  all  proportional  to 
An  or  At. 

II.  A  Cubic  Test  Body 

In  order  to  see  in  detail  and  more  easily  the  mechanism  causing  those 
quantities,  An,  <Ap>,  <  ApAp>,  etc.  ,  we  suppose  the  test  body  to  be  a  cube  with 
edge  length  (side  length)  D.  We  take  the  rectangular  coordinate  axes,  x,  y,  z,  as 
perpendicular  respectively  to  three  paris  of  faces.  By  assuming  that  the  direction 
of  v  is  perpendicular  to  one  pair  of  faces  which  is,  for  instance,  perpendicular 
to  the  x-axis 

v  =  (v,  0,  0),  (C-23) 


it  is  easy  to  calculate  those  quantities  of  our  present  interest.  The  face  which 

faces  to  the  x-direction  is  denoted  by  S+x  and  the  face  in  the  -x-direction  by 

S-x.  All  the  quantities  related  to  S+^  are  denoted  by  symbols  with  subscript 

+  x,  and  so  on.  The  number  of  molecules  which  collide  on  the  side  S,  in  At  is 

+x 

given  by 


+  oo  c  -  v 
x 


n+x  "  J  /  J  (V'Cx)fd5 


Jc  A  tDr 


I 


z 

00 


c  -  o 

X 


m  .3/2 

Cxn(  T^kT>  eXp 


„  L(C  -v)3  +  C  2  +  C?  Jld^CAtD® 
2kT  x  y  z  1 


]jd" 
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.  „„  .  m  ,  1  /  2 

'  At  a  -  'TSEt  > 


00 


m 


JC=0C»'1+  2M-  2  C*  V| 

L-  X 

m 


X  CXP  (  '  IkT  CxS  >  d  C, 


AtD*n<  JS^/2  jj<T_  +  _2L  (  2kT  3/2{ 

2rrk  T  (  m  4  kT  m  I 

|  2^tt  m 


2  *  f 


(C-24) 


where,  as  before,  mv  2/kT  <  <  1  is  assumed.  The  number  of  molecules  which 

collide  on  the  opposite  side,  S-x,  is  given  by 

I  1  .  2kT  ,1/2  v  I  ._ 

An  =  At  D2 nf  -  /=-  (  - )  -  —  >  (C-25) 

-x  I  2/n  m  2  \ 

We  see  that  the  number  of  particles  colliding  on  S+^  is  larger  than  those  on 

S-x.  On  the  sides  perpendicular  to  the  y-  and  z-axes  are  respectively 


An,  =An  =An,  =An 
+y  -y  +z  -z 


A.  _o  1  ,  2kT,l/2 

=  At  D-n  ( - 

2^n  m 


(C- 26) 


The  total  number  of  those  molecules  is 

3  .  2kT  ,1/2  A< 

An  =  -=  D  n  ( - )  At 

^tT  m 

(i)  <Ap  >  .  The  total  momentum  transferred  on  S  +  x  is 


(C-27) 


c  c  =  +  oo  c  =  v 

y«»  z  x 


<  Ap  > 
x  +  x 


-n 


j  2  m(v-  c^)p  fd3  c  At  D2 


c  ,  c  =  -  oo  C  =  -oo 
y  z  x 
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m  1/2 

=  At  D2  (-2m)  n  (  )  x 


oo 


J  <V  <  1  +  2  Cx  V)  X  CXP  ('  *2kT  CxS)  d  °x 

C  =  0 
x 


2mnD®  At 


J-  <^i  ♦»< 

4  m 


2kT  v1  /2  1 

m  J 


Similarly  we  have 


<Ap  >  =  + 

x  -x 


2mnDaAt  JtT  .  2kT 

^  4  m 


tt  ,  2kT  .  2kT  .4 
; —  ( — ~  )  -  v(-— > 
4mm 


Hence 


<Ap* 


> 

At 


<  Ap 


> 

x  4X 


<  Ap 

x 


> 

-X 


4mniy  At  .2kT.l  /  2 

jrr  '  m  ’ 


-  <Apy  V=  <Apy  %  =  ‘  <&p*  >+*  =  <Ap*  >*= 


kT 

mnDa  At  (  —  ) 
m 


and 


<Ap  >  =  <  Ap  >.  .  =  0 

Fy  At  rz  At 

Summarizing  these,  we  may  write 

-  4 

<Ap>4i  =  -  T  mv  An 
At  i 


(C-28) 


(C-29) 


(C-30) 


(C- 31) 


(C-32) 


(C-33) 
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We  see  the  cause  of  <Aps^  *n  the  unsymmetry  appearing  in  An+  and  An 


and/or  in  <Ap  >,  and  <Ap  > 

X  +X  X  -X 


(ii)  <Ap  Ap  >.  Since  the  collisions  are  binary,  we  may  write 


<  Ap  Ap  >  =  <Ap  Ap  >,  +  <Ap  Ap  > 

x  x  At  x  x  +x  x  x  -x 


<Ap  Ap  >, 

X  X  +X 


=  JIfx 


C  =00 
X 


4ms  (v-cx)3  fd3c  At  D2 


Similarly 


MDa4mSn<lSkT)1/2x 


x  exP  (-  C  3  )  d  C 

r  2kT  x  x 


J  v* 

C  =0 
x 


m  _  _  . 

2kT  2  °x  V) 


-  d>  » ( 2-a.T,1/2x  i  (-*£>%  £  v  5/2 


=  4m3  At  D?  n  — 
. 


1  2kT.  3/2  3  .  m.l 

+  T  V  (2kT> 


.  .  2  .  .  ,^a  1  .  2kT.3/2  3  .  m  .1 

<Ap  Ap  >  =  4m2  At  D  n  — =•  (  - )  -  —  v  (  ) 

*x  *x  -x  2  Jtt  m  4  2kT  | 

. 


Hence 


<Ap  Ap  >=  4m2At  D*n  (  2^^  )^2 
x  x  m 


4  a  /2kT.  . 

=  T  m  An 

5  m 


(C-34) 
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Similarly  we  have 


4  ?\tT 

<Ap  Ap  >  =  <Ap  Ap  >=  —  ma  ( - -  )  An 

x  y  x  z  3  m 


(C- 35) 


Since  Ap  and  Ap  are  due  to  two  groups  of  particles  which  are  independent  of 
x  y 


each  other,  we  have 


<Apx  Apy  >=  <Apx  >  <Apy>  =  ° 

(C-36) 

<Ap  Ap  >=  <Ap  ><Ap  >=  0 
y  z  y  z 


(iii)  <  Ap  Ap  Ap  >.  By  taking  advantage  of  easy  treatments  regarding 
this  simple  test  body,  we  may  have 


<Ap  Ap  Ap  >  =  <Ap  Ap  Ap  >  <Ap  Ap  Ap  > 

x  x  x  At  x  x  x  +x  x  x  x  -x 

<Apx  Apx  Apx  =  III  ^  2m(v‘  cx)  3  <v*  Cx)  fd3c  **  02 


[J 


■  ■  •**  D’8“>''<-rSkT,‘,2x 

00 

r  C  4  (  1+  ~  2  C  v)  exp  (-  ■—=  c  s)  dC 
Q  x  2kT  x  r  2kT  x  x 

X 


.  .  -.jn  3  1  3  i— .  2kT.5/ 2  .  m  .  _  ,2kT  .3 

=  -  A  t  DH8m3  n  gjTrf— )  +  (  ^  )  2  ™(  — ) 

=  -  At  D2  8m3  n  ^  3 


,  m  1/2 

X(W 


3  r-  2kT.  1  / 2 
a  ,»TT  <  ~)  +  2  v 

o  m 


.  2kT,3/2  1 

( - )  x 

m  Jn 


Similarly 


<Ap  Ap  Ap  >  =  +  At  D28m3n 

X  X  X  “  X 


3VtT  ,2kT.l/2  , 

8  -  2v 


,2kT.3/2  1 

x  /r 

m  JTT 
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Hence 


3.  ,  2kT\3/2 


<Ap  Ap  Ap  >  =  -  32  At  Da  m3n  ( 

rx  x  At  m 


r* 


32  3  .  2kT.  . 

—  m  (  - )  v  A  n 

3  m 


(C-37) 


Similar  calculations  show  that 


<  Ap  Ap  Ap  >  =  <Ap  Ap  Ap  >  =  0 

y  y  y  At  z  z  z  At 


(C- 38) 


All  the  other  components  where  Ap.  and  Ap^  (i  i  j)  appear  as  mixed,  vanish. 

The  reason  is  that  each  single  collision  is  independent  and  each  single 

collision  does  not  contribute  both  Ap.  and  Ap^  unles  i  =  j. 

(iv)  <Ap  Lp  Ap  Ap>.  The  method  of  calculation  is  similar  to  the 

preceding  cases.  It  is  simply  noted  that 

<Ap  Ap  A  ip  Ap  >.  =  "^Ap  Ap  Ap  Ap  > 

x  x  x  x  At  y  y  y  y  At 

=  <Ap  Ap  Ap  Ap  >  (C-39) 

Z  Z  Z  Z  "it 

oc  At 

The  other  components  vanish. 

Through  those  investigations,  it  is  obvious  that  a  molecule  with 
c  =-  |c  |,  relative  to  the  laboratory  has  velocity  component  ^x="lcxl  *v 
relative  to  the  test  body,  and  is  reflected  from  S+x  with  velocity  |c^|  +v 
relative  to  S+x  or  |c  |+2v  relative  to  the  laboratory.  On  the  other  hand,  a 
molecule  c  =+  |c  |  has  velocity  C  =  |c  |-v  relative  to  the  test  body  and  is 
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reflected  from  S-x  with  velocity  -  lcxl+v  relative  to  S-x,  or  -  |cx|+2v 
relative  to  the  laboratory.  We  see  that  the  test  body  accelerates  molecules 
towards  the  direction  of  v  and  is  accompanied  by  a  wake. 

We  may  also  note  that  the  number  of  molecules  colliding  on  S+x  in  At  is 


An,  =  fx  (c  -v)  dc  dc  dc 
+x  J "  x  x  y  z 


where  the  domains  of  integration  are 


c  -v:  -  oo  to  0, 
x 


:  -  oo  to  +  oo. 


On  the  other  hand 

An  =  fx  (c  -v)  dc'  dc  dc  . 
-x  JJJ  x  x  y  z 


c  -v  :  0  to  +  oo, 
x 


00  to  +  00 


Here  all  the  velocity  components  must  be  given  related  to  the  laboratory. 


56 


APPENDIX  D 


MARKOFF'S  METHOD  OF  RANDOM  FLIGHTS 


Markoff's  method  of  random  flights  is  summarized  for  the  convenience 
4 

of  application  in  this  article. 

An  elementary  momentum  given  to  a  test  body  is  p^  with  the  probability 
of  appearance  in  unit  time  w^,  Here  we  suppose  that  p^  and  w^  are  functions  of 
a  set  of  independent  variables  q}  ,  q2  ,  . . .  ,  .  Also  we  assume  that  there  are  v 

of  such  elementary  momenta. 

Let  us  denote  the  probability  of  the  total  momentum  being  between 
£  and  P  +  dP  by 


W(P)  d3P 


Then 


W(P)  d3P  =  J---J  A  n_  j  wr(q)d8q 


(D-l) 


Here  A  is  a  function  which  satisfies  the  following  conditions 
A  =  1  whenever 

v  _ 

p-  5  dp<l 

r  =  1 

A  =  0  otherwise. 

It  is  known  that  A  is  given  oy  Dir  ichlet's  integral 

,  1  fT  sin  (ydpx  Px)  sin  (  j  dPy  Py)  sin  (  \  d  Pz  p,) 

A  =  n—  JjJ  - i -  - 2 - 


(D-2) 


=  JJJ 


x  exp  L  i  P  .  (  v  p  -  P)  jdp  dp  dp 
>  r  x  y  z 

r=  1 

(fhp  CXP  [iP'  {l  Pr  *  P)  jd"P 


(D-3) 
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Hence 


W(£)  =  ~j3  J---J  exp  (-ip.  P )  A  (p)  d3p 

A(p)  =  ft.  f  Cexp  (ip.  p  )  ]  w  d8q  (D-4) 

r- i  J  r  r 


If  p  and  w  , 
r  r 


as  functions  of  the  q's,  do  not  depend  on  subscript  r. 


we  may  have 


A(p)  =  |  J  C  exp  (ip  .  p)  ]w  d®q  |V 
For  application  in  case  of  V  >>  1,  we  may  write 

A(p)  =  j  1  -  —  ^  v  *v  j*  ^  exp  ^  ip  •  P)^  w  d8<l| 

j  v  -  v  J  L  exp  (ip  .  p)  ]  w  dsq  | 


(D-5) 


=  exp 


(D-6) 


APPENDIX  E 

MULTIPLE  COLLISIONS  OF  AN  ELASTIC  TEST  BODY 
IN  A  RAREFIED  GAS 

As  shown  in  Section  II  of  Part  II,  if  is  much  larger  than  unity, 


we  have  to  calculate 


VAt 


1P>4.=  i.  %Pr 


r  =  1 


=  l-~l  ,rj  ____r  . 


(II.  2.  26) 

Since  <  A p  >  is  linear  with  respect  to  p  ,  there  is  no  difference  in  <Ap  >. 

At  rx  At 

between  the  case  of  binary  interaction  (v  <1)  and  the  case  of  multiple  interaction 


For  calculating  <Ap  Ap  >t>  in  the  approximation  of  neglecting  the 
terms  of  0  [  ^  ],  we  may  assume  that  the  velocity  of  the  test  body  relative  to 

the  gas  is  zero,  as  shown  in  Appendix  C. 

Case  Il-a.  The  test  body  is  a  spherical  body  submerged  in  a  gas  in 
thermal  equilibrium. 


On  an  elementary  surface  in  Fig.  5, 
dS  =  sin  8  d  8  dcp  R^, 


we  consider  a  local  rectangular  coordinate  system  (§,ti,C)  where  the  C-axis  is 
inwardly  perpendicular  to  dS  and  the  ?-and  ri-axes  are  parallel  to  the  surface. 
The  distribution  of  the-field  molecules  is  given  by 

f  =  "  (  2TkT,3/2  CXP  [  *2kT  (C-E  +  S?  +  CQ  ]  3  {E'l) 

The  number  of  field  particles  which  collide  on  dS  in  At  is  given  by 
dw  =  fc^d3.cdS  At, 

where 

d  S  =  sin  0  d  8  d  c  Ra, 

^  ;  0  to  2  tt,  (E-2) 

8  ;  0  to  n . 
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One  may  obtain  <v>^t  directly  from  (E-2) 

r 

<  v  >  =  I  d  w 

At  J 


=  4  ttR3  n 


(E-3) 


A  particle  gives  the  test  body  momentum  p 


p  =  2  c  .  mR/R 
r  C 

The  total  momentum  is 
<V> 

-  c  At  - 

P\l,  p- 


(E-4) 


(E-5) 


According  to  Markoff's  method  given  in  Appendix  D,  the  probability  of  P 


being  between  P  and  P  +  d?*  is 


1 

W(P)  d^P  =  ^3  JJJ  exp  {-ip  .  P)  A  (p)d5p  d3  P 

-  OO 

A(^  =  [  J  exp  (ip.  p  )  ]  V  ^.t 

At 

Since  A(p)  is  expected  to  be  independent  of  the  direction  of  p  (spherical 


(E-6) 

(E-7) 


symmetry)  we  first  calculate  A(p)  by  taking  pin  the  direction  of  the  z-axis. 


A(d)  =  [  JJJ  exp  (i  p  .  2m  c-  cos  9  )  - ^ 


fd3c  c-  RH  sin  9  d  9  d  cpAt 


<v  > 


At 


<  v  > 


At 

.  [  4nR^n(m/2nkT)i/ZAt  J  8in(^cg9) 
<V>At 


2mp 


,  m  _ .  .  T  <  v> 

exp(-TTrcC  ',dcC3  At 


By  expanding  sin(Zmc^  p)  in  Taylor  series,  we  have 


oo 


j 


sin(2m  c^  p) 
2mp 


exp(-  CV  dcc 


kT  4  1 A 

— -  L  1  -  7  mkT  p2  +  77  (mkT)"  P4 
m  3  15 


2 


3x5x7 


(mkTf  p'  +- 
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Considering  that  <v  ^  >>1,  we  obtain  for  A(p) 

A((T>=  exp  [<v>At  (  -  |  mkT  ps  +  j-|  (mkT)sp4  -  - )] 

Remembering  that  A(p)  is  independent  of  the  direction  of  p,  we  write  for  A(f^ 
A(p)  =  exp  [  -<  v  >  |  mkT  fP  ]  x 

x  j  1  +  ||  <v>t  (mkT)s  p4  +  ---  j  (E-8) 


By  remembering 
oo 


exp  (-  a3xp)  .  sin  fix  .  x  dx 


_  ..Vzi..  exT}  / .  _S!_  , 

-  4  a3  exp  (  4a?  ' 

r°  .  Oa4B-  20cf  fi3  +  B6) 

exp  (-  ax")  •  sin  fix  .  xB  dx  =  - K—L 


6  4  a11 


we  finally  obtain 

W(P)  =  “ ^  HI  exp  (-i  d  P  cos  6)  A(d)  p2  dr  sin  fid  0  d  r 


X  exp(_4a*") 


i  00 


P  F 


8n^(<v>(  f  mkT) 

At  J 


7/2  +  0  Hr—  )7/Z] 


<  M> 


fit 


3PS  1 

XCXP(-  1  6 ' <  v  >~ mkT" 


n 


We  may  ignore  the  higher  order  terms  of  (  — ).  It  is  shown  that 


W(P)  d  P  =  1 


(E-  10) 
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We  calculate  first  <  Ap  Ap  >  : 

At 


<  Ap  Ap  >  +  <  Apv  Ap  >  +  <  Ap  Ap  > 
xx  y  y  z  z 

=  jyj ps  d3p 


(E-ll) 


=  8  mkT  <v>  , 

At 

Since  W(P)is  independent  of  the  direction  of  P,  it  is  a  simple  matter  to  show 
that 


<  Ap  Ap  >  =  <Ap  Ap  >  =  <Ap  Ap  > 
x  x  y  y  z  z 

=  f  <V>  t  mkT 


(E- 12) 


<Ap  Ap  Ap  >  :  It  is  conceivable  that  <Ap  Ap  Ap  >is  proportional  to  v. 
It  is  not  a  simple  matter  to  calculate  <Ap  Ap  Ap  >  according  to  the  Markoff's 
method,  since  A  (p  )  depends  on  the  direction  of  p,  We  write 


instead  of  (E.  2), and 


dw  =  f  (c^  +  v  cos  0  )  d3c  ds  At. 


p.  =  2  (C  -  +  v  cos  0  )  m  R/R 


instead  of  (E.  4)  ,  by  taking  v  in  the  direction  of  the  z~axis. 


P  ■  p.  =  o  IpJ  cos  S 

where  iis  a  function  of  £,  -.  and  9  ,  -  ,  and  the  manipulation  must  be  complicated. 

D  r 

It  is  a  simple  matter,  however,  to  calculate  «.*  'p^  Ap  Ap  Ap  >,  by  ignoring  the 
terms  of  0  Lv9]by  means  of  W(f^  given  by  (E.  6) 


<  Ap  Ap  Ap  Ap  >=  <Ap  Ap  A p  Ap  > 

X  X  X  X  ry  ry  *y  ry 


=  <Ap  Ap  Ap  Ap  > 

z  z  z  z 


(P  cos  9) 4  W(P)  d3P 
r  oo 

=  2  tt  j  cos4ede^P6  W(P)  dP 


It  is  easily  shown  that 
oo 

f  P6  W(P)  dP  oc 


<  v>, 


7/2 


At 


=  <v>; 


<v> 


3/2 


At 


At 


(E-13) 


Cubic  Test  Body.  If  the  test  body  is  cubic,  and  one  of  the  faces  is 
perpendicular  to  the  direction  of  v,  the  calculation  is  very  simple.  Here  An  and 
<AP  >  do  not  change  from  those  obtained  in  the  case  of  binary  interactions 


A  n  =  AtD^n(-^I)1/2 

/tT  m 


<Ap'>--  —  mv  An 


(E-14) 


Case  II- bj.  The  particles,  the  number  of  which  is  given  by 


v  =  J|  fd''  c  c^  tx  6D2, 

are  distributed  at  random  pn  the  6  faces.  Because  of  the  steady  motion  of  the 
test  body  perpendicular  to  one  side,  such  as  shown  in  the  figure, 


Fig.  6 
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the  number  of  molecules  which  collide  on  S,  is 

+x 


+QD 


C  =  V 
X 


.U  f  DMx(-Cx+v)d3c=  An+x 


^  -  00  c  =  -  oo 

X 

z  v  >  0 


and  the  number  of  those  which  collide  on  the  opposite  side  is 
+oo 


jj 


CX  =  +  oo 

f  D2  f  x  (c  -  v)  d  3  c  =  An 

0  X  -X 

c 

X  =  V 


=  -  00 


By  ignoring  the  members  of  the  order  of  9  We  have 

+  oo  +oo  o 

Vx=J  I  I  D?fx(-cx+v)d3c 

c  =  -  oo  c  =-oo  t  =  -  CD 


An, 


=  Dan  f  (*c-„  +  «P  (-  2kT  V*  dC- 


2  nkT 


m 


i 


C  =  -  00 
X 


=  D?n  (  ) 

'  2  ’"kT 


2kT  +  v  VtT  (  2kT  j 
2  m  2  m 


At 


1  k  T  F 

=  D'"W/2(-sr» 


*  T  1 


An  =  D  ~  n 
-x 


(±T  v 

m 


v 

2 


At 


Similarly, 


An  -An  =  '  n  =  f n 

+y  -y  +z  -z  (2n)S  m 


D~n  ,  kT.5 
S  (  — )  At 


and 


„  M  1  ,  kT  J. 

An  =  6  D3  n  - — r~i  (  -  )  't. 

(2  n)  £  m 


(E- 15) 


(A 


The  momentum  due  to  An+x  in  the  negative  direction  of  the  x-axia  is 


x  =  v 


Ap 


-  v)  d3C 


*x  ■  J'J|  J>‘  ‘  <-  cx  +  v)  2m  |c> 

X 

=  nDMI?iT,42mJ0  ^  +  2  Cx  eX^  (*  2kT"  Cx>,d 


C  =  -  CD 
X 


n D2  {  -r;J?x2m 
1  2  rrkT 


^rT  ,  2kT.3/2  ,  2kT 

—  ( “5T>  -  2v 


2m  Vn”  ,  2kT  ,  ,  2kT  .j 

n  D2  ■=  -  — —  (  —— )  -  v  ( ) 

tt  4  m  m 


Similarly 


Ap 


+oo 

•**15. 


D2f  (c  -  v)?  2m  d3  c 


=  n  D°  <  + 

<r 


2m  i.  /T  ,  2kT .  , 2kT. s 


(  -  v  (  — ) 

m  m 


Hence 


AP  =  (Ap+  x  +  AP_X  ) 


V 

V 


nD2  4m  ,  2kT  *• 
— jF==~  (  — — -  )  v 
Jtt  m 


—  m  An  v 


(E- 16) 


By  taking  the  Maxwell  distribution  for  f  we  have 

3  ..2  ,  2kT.s  — 

v  =  -==  D2n  - )  T 

JrT  m 


(E-  17) 
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As  before,  we  calculate  A(p)  obtaining 

2mc, 

A(p) 


exp  L  ip  — ~  ]  +  exp  L  -  i  p 


2mc , 


X  T 


2mc  2mc. 

+  exp  [  +  ip  — y— ]  +  exp  L  -  ip  — -  ] 


2m  c 


c 


2  me' 


+  exp  [  +  ip^  — - -  ]  +  exp  [  -  ip 


L 


11 

It 


f d3c  c  n8  T 


•  1  .  2m  ?  3 

2  '  —  '  V 


p8  +  p8  +  p8  "I  fd3c  c -6  D8  t 


2  P8  c,  8 


!  .  i  (  2-2-, 

2  '  t  '  3 


fd3c  c,  6  D8  T 


W(F)  = 


4  P  8 

exp  -  -  vmkT  (  — ) 


1  rr;'  f 

JJJexp('iFxPx)exp  l- 


t  » 


y  mkT 


r  V  4  1 

x  exp  (-iFy  py)  exp  ^  v  3  mkTJ 

r  p22  4  1 

x  exp  (-iF  p  )  exp  I-  v  — 5—  *r  mkT 
z  z  L  t  3  J 


xdp  dp  dp 
x  y  z 


1 

" 

F8 

’  V 

U< 

2m  kT' 
t  rn 

3/2  exp 

* 

i  v(2-fr^X 

jJJ  V  W<F>  d3F  = 


+  OO 


_f  V  2m  3  kT  i  J  Fx‘  exP 

2JTI  y  (  — )  (—  )]  .  a, 

2  v  2m  ■  k  T 

3  '  1  m 


4  v  .  2m.  ?  kT 
T'  T  '  m 


(E.18) 


—  dF 

J  X 

(E- 1 9) 
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Hence 


8  8 

<Ap  Ap  >=  —  v  mkTAt/T=  —  mkT  An 
x  x  3  3 


( E-20) 


Here  again  we  see  that 


<.ip  Ap  Ap  Ap  >oc  v! 

X  X  X  X 


(E-21) 


Case  II- bs .  On  each  side,  the  number  of  particles  is  assigned  as 

v  -»  -* 

—  due  to  some  physical  condition.  In  this  case,  we  have  to  calculate  <Ap  Ap  >on 

each  side.  For  example,  by  taking  the  side  which  is  perpendicular  to  the  x-axis 

on  the  positive  side,  we  have  only  one  finite  component,  <Ap^  Ap  .  For 

A(p)  we  have 


A<Px) 


2m  kT 

m  '  x 


„  v  __  m  2kT.£ 


1  v  2n 
+  =  exp  j  -  T  (— 

w,Fx>  =  -hr  N-*r.v  A<v d °x 


Noting  that 


1 


ItT  v  Zm.g  kT  T" 
6  t  '  m 


exp 


fF, + 

4  v  2m  3  kT 
6  (  t  '  m 


(E-22) 


“  t  '  m '  J  TT 

— r-r - 5 — r •  =  —  v  >>1, 


(E-23) 


.  ,2m*’  kT 

4 v  (— )  — — 

t  m 

F  1  =  F  +  v^rT  —  ( - )5  ranges  from  -do  to  +  oo  as  F  changes  from  -oo  to  0. 

X  X  T  X 
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Hence  we  have 


9° 


i  F.  W  (F  )  dF  =  -  JT  —  ~  ( 
J+x  x+x  6t 


v_  rn  ,  2  k  T.£ 
m  ' 


f  F  3  W (F  )  dF  =  ^  (2j2-)  —  +  2tt  (-£-)*(  Zf-)1  — 
J  +x  x  +■  x  3  T'm  6  T  m 


{E  24) 


and  for  the  opposite  side 

r  F  W(F  )  dF  =  /tT  —  (  —■)  (— )^ 
J  -x  x  -  x  *  6  t  '  1  m 


oo 

{•  r> 


i  F  W(F  )  dF  =  ~  <r~) 
J-x  x-  x  3  t  rn 


v  ,2m\  2kT  +  2  *  {—)  —  (E  25) 

6  t  m 


For  calculating  the  total  in  the  x-direction,  we  put 


A  t 


<Ap  Ap  >  -  TS<(FX  +  F  )3  >x  — 

x  x  At  tx  -x  T 


=  T  At  (  <  F  ,  >+<F  >  +  2  <  F ,  F  >) 

+X  -X  +X  -X 


Since  F  and  F  are  independent  events,  we  nave 
+x  -x 

<  Ap  Ap  >  -  t  At  (  <  F  ~  >  +  <  F  3  >  +  2  <  F  ><F  >) 

x  x  At  +x  -x  +x  -x 


2  ,2m  3  kT  At 

=  3  v(7>  -77T 


-  mkT  An 


(E  26) 


Case  II- b3  .  We  now  assume  that  the  distribution  of  field  particles 
in  the  configuration  space  is  at  random.  It  is  assumed,  however,  that  each 
particle  has  the  same  momentum  components 
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lpix1=  'Piy1=  'Pizi  = 


m  |  c  I  fd3c 
x 


!«■ 


nn  ,  m  t  .  kT  . 
-  m(2~T)  (~) 


1  - kT  .  a 

(  2  tt)  |  m  (  m  ’ 


mkT  £ 
'  2n  ' 


(=■27) 


In  this  case  we  obtain 


<v>t=  6nD^t(^)^/m 


c iZL,  * At 

(2tt)  g  '  m  ’ 


(E  28) 


<ap  >  =  Drn  ~  f  +  v 

x  At  ) (2 tt )  l  m 


At 


,  ,  mkT ,  5  .  . 

x  {  (  y^~)  +  m  v  I  x2 


,  mkT  ,g  , 

(— )  -  m  v  }  x  2 


-4 


1  ,kT  l  ,  .  mkT,?  .  _  „ 

m  +  *  vMD  n 


8(f  v  A t  D2  n 


3  m  <V>At  V 


(E  29) 


In  this  case 


Case 

signed 


A(P  )  =  <  2  cos 


x  '  2 


TT 


2  P  (~rn+2co. 


[2“z  'TT"1']  ( 


1  V  ' 


2  p  <*£>’ 
y  2n 


,  ,  4mkT.  1  ,  P2 . 

=  exp{  -  ( -2"^ — )  2  V 


j  a  mkT  } 

exp  J  .  vp*  —  j 


1  ?  00 

(PX)=  2T  J  exp  (-i  Px  Px)  exp 

-CO 


V  p 

X 


mkT 

3tt 


d  p 


(E-30) 


+oo 


=  2n  J  cos  (Px  p  )  exp 

_O0 


mkT  n  9 

-  v  TT"  px* 
3  TT 


dp_ 


2  TT 


/TT 


(v 

'  3  TT  ' 


exp 


.  P  8  1 

x 

"  mkT 

L<v—  J 


i 

2  TT' 


(  V 


[-  PXR 

77^T 

3  n 


(E-31) 


<ApXAPx>  =  iV* 


,  mkT 
<v  — 1 


£1  ,4  v  .2*1) 

2  '  3tt 


3  TT 


v  mkT 


v  =  <  v  > 
At 


(E-  32) 


II-b4.  Each  particle  exerts  the  average  force  and  —  particles  are  as- 
on  each  side.  It  is  obvious  that  Ap^  has  no  fluctuation  and  <Ap  Ap>  vanishes. 
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APPENDIX  F 


as 


USEFUL  INTEGRALS 


The  integrals  which  often  appear  in  the  present  manipulation  are  listed 
follows  7  ,  8  i; 


oo 


j  X2n  exp  (  -X  X2)  dX  =  Hhr'l0  ^1*  <  "IZ+1 - ^  ’ 


n  ! 


I  X2n+1  exp  (-X  Xs)  d  X  = 

J  _  .  n+1 

o  2  \ 

oo  i 

j  exp  (-XX2)  dX  =  |  (-jr  , 


00 

r 


J  -  exp  (-X  X?)  dX  =  -J  (p-)*  , 


j  X  exp  (-  XX^)  dX  =  — 

o  * 

fx» 

V 

O 

<3°  •  1 

j  X3  exp  (-IXs)  dX  =  jp- 

o  * 

OO  J 

J  X4  exp(-XX2)dX  =  |(p)r 

o 

00  1 

f  Xs  exp  (-XXE)  dX  i  -p- 

o 

oo  cl 

j  Xs  exp  (-XX~)  dX  =  — 
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Jx7  exp  (-xx2)  dX  =  pr. 
o 


5  2/ 

Jexp  <-p2  X2)  cos  (qX)  dX  =  y  “  exp  [  -q  A  p  2  ] 
o 

00  ^2/ 

Jexp  (-p2  X2)  sin  (qX).XdX  =  exp  [  -q  f*  P2] 

o  ' 

Jexp  (-p2X2)  cos  (qX).  X2  dX  =  ~  exP  [  ^  p3] 

o 


exp  ,-p»x=,  .in  <qX>.  X-  dX  ^ 

f  exp  (-P»  x=|  CO.  (qX).  X*dX  .  f^yvaa-g*^  [-q^p=] 


j 

o 


"exp  «-p<  X=)  .in  «*..  X»dX  -  exp  t  -  ,^3 
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